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ABSTRACT 

Pressure  pulses  caused  by  irradiation  of  a model  gallium 
arsenide  solar  cell  with  a Carbon  Dioxide  TEA  Laser 

were  examined  using  power  densities  of  the  order  of 

7 2 

10  watts/cm  . The  pressure  pulses  were  monitored  with  a 
Sandia  type  quartz  pressure  gauge.  It  was  discovered 
that  the  relatively  low  power  densities  used  were  capable 
of  removing  the  silicon  dioxide  antireflective  and  gold 
contact  layers  of  the  solar  cell  after  only  a few  shots 
of  the  laser.  An  exponential  relationship  between  the 
initial  thicknesses  of  these  layers  and  the  pressure 
pulse  generated  in  the  gallium  arsenide  substrate  was 
indicated  for  gold  contact  layers  of  less  than  5000  A 
thickness.  Evidence  was  found  that  the  principal  pressure 
generation  mechanism  is  thermomechanical.  Gold  films  of 
thickness  greater  than  5000  A were  found  to  be  able  to 
absorb  the  power  densities  used  with  no  apparent  damage. 
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I.  INTRODUCTION 


The  development  of  solar  cells  as  versatile  and  dependa- 


ble sources  of  power  for  satellites  is  well-known.  Virtually 


every  long-range  space  probe  launched  by  the  United  States 


or  the  Soviet  Union  in  the  last  decade  has  been  powered  by 


photovoltaic  devices.  Solar  cells  provided  the  power  that 


operated  the  television  cameras  of  the  Ranger,  Surveyor  and 


Lunar  Orbiter  spacecraft  that  made  possible  detailed  photo- 


graphic studies  of  the  surface  of  the  moon.  These  devices 


provided  the  power  for  the  Mariner  probes  of  the  planet 


Venus.  The  recent  Skylab  vehicles  also  utilized  solar  cells 


in  their  electrical  generation  systems.  The  application  of 


solar  cell  technology  has  led  to  the  development  of  the 


systems  of  communications,  navigational  and  meteorological 


satellites  currently  being  employed  by  the  United  States 
and  other  countries  [Ref.  1], 


The  general  construction  and  characteristics  of  solar 


cells  are  treated  in  detail  by  several  authors.  A good 


review  of  these  topics  is  contained  in  Direct  Energy 


Conversion  by  S.  W.  Angrist  [Ref.  2].  Of  particular  note  in 


this  treatment  is  the  requirement  that  for  maximum  theoretical 


conversion  efficiency  the  solar  cell  material  should  have  an 


energy  gap  of  between  1.1  and  2.3  electron  volts.  Gallium 


Arsenide  (GaAs)  with  an  energy  gap  of  1.35  ev.  satisfies 


this  condition. 
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A significant  amount  of  research  has  been  undertaken 
previously  to  ascertain  the  effects  of  laser  irradiation 
on  solar  cells  and  phototransistors.  A study  of  irradiation 
of  phosphorus -doped  silicon  solar  cells  by  a pulsed  laser 
! performed  by  Giuliani  [Ref.  3]  identified  two  permanent 

2 

damage  effects:  Below  laser  energy  densities  of  5 joules/cm 

no  permanent  damage  to  the  phototransistor  was  observed. 
Irradiation  at  these  low  energy  densities  was  found  to 
increase  the  photoconductive  response  of  the  device,  however. 
At  energy  densities  above  10  joules/cm  permanent  changes  in 
the  electrical  characteristics  in  the  form  of  increased 
dark  current  were  observed.  Examination  of  phototransistors 
damaged  by  laser  energy  densities  in  this  range  enabled 
'Marquardt,  Giuliani  and  Fraser  [Ref.  4]  to  explain  this 
increased  dark  current.  Evidence  was  found  that  the  melting 
point  of  the  silicon  had  been  exceeded.  A model  for  phos- 
phorus ion  migration  was  developed  by  treating  the  phosphorus 
as  a dilute  impurity  in  a melted  silicon  matrix  and  then 
analyzing  the  migration  of  the  phosphorus  ions  as  a diffusion 
problem.  The  subsequent  refreezing  of  the  silicon  then  traps 
the  ions  in  their  new  distribution  and  causes  the  increased 
dark  current  observed. 

As  evidenced  by  the  above  work  a mechanism  explaining 
the  changes  in  the  electrical  properties  of  phototransistors 
under  laser  irradiation  has  been  developed.  However,  the 
mechanical  damage  done  to  solar  cells  by  laser  radiation 
has  not  been  so  thoroughly  examined.  By  investigation  of 
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the  momentum  delivered  to  a GaAs  solar  cell  by  laser  irradia- 
tion it  was  the  intent  of  the  current  investigation  to  amplify 
the  knowledge  in  this  area.  A model  of  a GaAs  solar  cell 
was  constructed  as  shown  in  Figure  1.  The  thickness  of  the 
gold  film  was  varied  and  a measurement  of  the  transmittance 
of  the  gold  film  was  made  at  the  wavelength  of  the  Carbon 
Dioxide  laser  used  (10.6  urn).  This  was  done  to  determine 
whether  the  gold  coating  would  absorb  a significant  amount 
of  the  incident  laser  energy.  The  solar  cell  model  was 
then  constructed  and  irradiated  with  light  from  a CC>2  TEA 
laser.  The  impulse  transmitted  to  the  GaAs  was  measured  as 
a function  of  time  and  of  the  thicknesses  of  the  gold  and 
SiC^  film  using  a Sandia-type  quartz  gauge  [Refs.  5 and  6], 
Finally,  an  experiment  was  performed  to  determine  the  energy 
threshold  for  removal  of  the  quartz  film  by  the  incident 
laser  beam. 


II.  PREVIOUS  WORK 


The  absorption  of  laser  radiation  by  a great  variety  of 
materials  under  a large  range  of  conditions  has  been  examined 
by  several  experimenters.  The  type  and  mechanism  of  damage 
is  dependent  upon  whether  the  target  material  is  opaque  or 
transparent  to  the  incident  laser  light. 

Damage  to  transparent  materials  by  laser  radiation  has 
been  the  subject  of  much  experimental  and  theoretical  work. 
Observations  of  the  damage  itself  are  plentiful  but  a clear 
understanding  of  the  mechanism  of  damage  has  yet  to  be 
developed. 

The  main  mechanisms  suggested  for  the  laser-induced 
damage  to  transparent  materials  are  [Ref.  7]: 

a.  Production  of  phonons  (hypersound)  in  a stimulated 
Brillouin  scattering  process  with  the  phonons  being  produced 
in  an  intensity  sufficient  to  lead  to  a fracture  of  the 
material . 

b.  Absorption  of  light  by  defects  in  the  material  with 

resultant  thermal  shock  and  microplasma  production.  Using 

a ruby  laser  with  20  psec  to  20  nsec  pulse  durations  and 

2 

depositing  laser  energy  densities  of  from  3.5  joules/cm  to 
16.0  joules/cm  on  the  target,  Milam,  et  al.  [Ref.  8]  demon- 
strated that  heating  of  inclusions  ranging  in  size  from  0.1 
to  5.0  ym  is  the  limiting  mechanism  for  damage  in  TiO^  and 
ZrO^  films  deposited  on  quartz  substrates.  A relationship 
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between  the  length  of  the  laser  pulse  and  the  size  of  the 
inclusion  heated  was  also  developed.  Short  laser  pulses 
can  more  easily  heat  small  inclusions  and  long  pulses 
preferentially  heat  large  inclusions.  It  was  shown  by 
Hopper  and  Uhlmann  [Ref.  9]  that,  in  the  first  approxima- 
tion, the  inclusion  that  is  most  readily  heated  by  a laser 

1/2 

pulse  of  duration  x^  has  a radius  proportional  to  (x  ) . 

c.  The  electron  avalanche  phenomenon.  This  phenomenon 
occurs  at  the  surface  of  a transparent  dielectric  [Ref.  10]. 
’'Free”  electrons  at  the  surface  of  the  dielectric  are 
accelerated  by  the  electric  field  of  the  laser  beam  sufficient 
ly  to  produce  ionizing  collisions  with  the  atoms  in  the 
surface  layer  thereby  producing  more  free  electrons.  When 
the  lattice  of  the  solid  acquires  energy  from  the  free 
electrons  so  produced  its  temperature  rises.  The  surface 
matter  then  melts,  boils  and  evaporates,  forming  a partially 
ionized  high-density  vapor  layer  above  the  surface.  This 
vapor  layer  absorbs  laser  light  strongly  and  the  absorbed 
light  energy  continues  the  ionization  process.  The  mechani- 
cal damage  resultant  at  the  surface  is  primarily  due  to  the 
plasma  formed  as  a result  of  an  optical  discharge  (thermal 
ionization  of  all  the  evaporated  matter)  and  a shock  wave 
generated  by  this  discharge. 

Damage  produced  in  transparent  materials  by  high-energy 
laser  beams  can  range  over  a wide  field  of  effects  extending 
from  microcracks  or  melted  voids  in  the  interior  of  the 


material  to  craters  and  cracks  on  the  surface  of  the  material 
[Ref.  11]. 

The  damage  effects  in  transparent  materials  are  highly 
power-dependent.  If  the  incident  laser  radiation  is  below 
a certain  power  threshold,  the  laser  beam  passes  through  the 
material  with  no  effect.  However,  if  this  threshold  is 
exceeded,  the  laser  beam  is  absorbed  and  damage  of  one  form 
or  another  takes  place  [Refs.  12,  13  and  14].  For  sodium 
chloride  crystals  irradiated  with  a Q-switched  ruby  laser 
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this  threshold  occurs  at  2.0  GW/cm  ; for  Silicon  Dioxide 
glass  the  threshold  is  at  470  GW/cm2  [Ref.  12].  Wang  [Ref.  14] 
has  measured  the  damage  threshold  for  copper  reflectors 
irradiated  with  10.6  ym  radiation  from  a CC>2  TEA  laser  to  be 
150  MW/ cm  . Experiments  [Ref.  11]  have  shown  that  the 
threshold  for  internal  damage  occurs  at  lower  levels  of 
irradiation  than  surface  damage. 

The  absorption  of  laser  energy  by  substances  opaque  to 
the  wavelength  of  the  laser  used  (absorbers)  has  also  been 
the  subject  of  considerable  investigation.  The  variance  of 
energy  absorption  with  irradiance  and  duration  of  exposure 
shows  three  distinct  regimes  differing  in  the  phenomena 
induced  in  the  material  by  the  laser  beam. 

The  first  of  these  regions  is  heating  of  the  material 
without  phase  change.  This  results  in  a temperature  rise 
in  the  material.  This  temperature  rise  may  lead  to  simple 
thermal  expansion  of  the  material  if  the  deposition  of  energy 
is  slow  enough.  However,  if  the  laser  energy  is  deposited 

14 


rapidly  enough,  say  in  a time  t comparable  to  or  less  than 
the  time  required  for  an  elastic  wave  of  velocity  v to 
cross  the  area  irradiated,  a pressure  wave  can  be  generated. 
That  is,  if 

t < a/v  (1) 

where  a is  the  radius  of  the  area  irradiated.  In  this  case 

the  irradiation  takes  place  so  quickly  that  the  material 

remains  "inertially  clamped"  during  the  absorption  process. 

Hence,  the  volume  of  the  irradiated  portion  of  the  material 

remains  constant  over  the  time  immediately  following  the 

- 8 

laser  irradiation.  A laser  pulse  of  duration  10  sec. 
satisfies  this  condition  [Ref.  15]. 

The  next  region  considered  is  laser-induced  melting  of 
the  target  without  vaporization.  In  this  region  the  damage 
caused  in  the  target  consists  of  melting  of  target  material 
to  a depth  dependent  upon  the  power  density  of  the  laser 
beam  irradiating  the  target  and  upon  physical  properties, 
such  as  the  heat  of  fusion,  of  the  target  material.  The 
optimum  pulse  length  required  for  this  effect  ranges  from 
50  ysec  to  5 msec  and  is  much  longer  than  the  pulse  lengths 
in  the  nanosecond  range  that  were  produced  with  the  lasers 
used  in  the  current  investigation  [Ref.  16]. 


The  final  region  in  which  absorption  by  opaque  materials 
takes  place  in  laser-induced  vaporization.  For  the  purposes 
of  the  experiment  performed  this  was  the  most  sicnificant 


region.  The  strong  energy  fluxes  concentrated  over  the  small 
area  of  the  focal  spot  of  a laser  beam  produce  such  intense 
evaporation  that  significant  pressures  can  be  generated  in 
the  absorbing  material  by  the  recoil  of  matter  ejected  by 
evaporation  [Refs.  17  and  18].  When  the  irradiation  takes 
place  in  a vacuum  a hot,  dense  plasma  is  formed  above  the 
material  and  the  recoil  of  the  particles  in  the  plasma 
transfers  momentum  to  the  target.  In  an  experiment  in  which 
aluminum  targets  were  irradiated  with  a Neodymium- glass 
laser  having  an  output  between  4 and  11.5  joules,  Krehl, 

Cooper  and  Schwirzke  [Ref.  19]  demonstrated  that  of  the  total 
mass  ejected  by  the  laser  beam  from  the  target  about  ten 
per  cent  was  ionized  and  that  this  small  ionized  component 
was  responsible  for  almost  80%  of  the  momentum  transferred 
to  the  target.  A simple  calculation  by  Askar'yan  and 
Moroz  [Ref.  17]  shows  that  the  recoil  pressure  of  the  ejected 
matter  is  many  times  the  radiation  pressure  of  the  laser  beam. 

Pirri  [Ref.  20]  has  developed  a comprehensive  theory  for 
the  momentum  transfer  to  a surface  from  a laser-supported 
detonation  wave  initiated  in  the  air  above  the  target  surface 
by  a laser  beam.  Although  this  theory  is  not  strictly 
applicable  to  the  vacuum  portion  of  the  current  investigation, 
it  is  described  here,  in  some  detail  as  an  excellent  descrip- 
tion of  the  mechanism  of  momentum  transfer  to  a surface  by 
a laser  beam.  The  model  is  based  upon  the  physical  processes 
which  occur  when  a high-power  axisymmetric  laser  beam 
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interacts  with  a surface  in  the  presence  of  an  atmosphere. 
Figure  2 illustrates  this  model. 

If  tp  is  the  laser  pulse  length,  at  some  time  t <<  rp, 
on  the  order  of  50  nsec,  a laser-supported  detonation  wave 
is  initiated  above  the  target  surface  by  breakdown  of  the 
small  amount  of  surface  material  vaporized  during  this  time. 
Once  the  detonation  wave  forms,  the  vaporization  stops  and 
the  force  delivered  subsequently  is  a result  of  shocking  the 
air  above  the  surface  to  a high  pressure. 

This  high  pressure  in  the  laser-shocked  air  means  that 
there  is  a significant  impulse  delivered  to  the  target 
after  the  laser  pulse  ends  as  the  pressure  relaxes  to  its 
ambient  value.  Hence,  the  dominant  time  scale  for  the 
momentum  transfer  process  is  tQ,  the  time  over  which  the 
impulse  is  delivered  to  the  target  surface  rather  than  just 
the  laser  pulse  length  Tp. 

Initially,  the  flow  of  the  hot,  dense  plasma  from  the 
surface  is  one-dimensional.  At  some  later  time,  t 2d,  the 
flow  becomes  two-dimensional  due  to  the  thermal  expansion 
of  the  plasma.  This  two-dimensional  time  is  of  the  order 


of  the  time  it  takes  the  leading  edge  of  the  plasma  to 
travel  a laser  snot  diameter  TRef . 211 . This  can  be  defined 


where  Dg  is  the  diameter  of  the  laser  spot  and  V is  the 

velocity  of  the  leading  edge  of  the  plasma. 

For  the  entire  pressure- time  history  of  the  momentum 

transfer  to  be  described  by  one -dimensional  analysis  the 

total  time  over  which  the  impulse  is  delivered  must  be  less 

than  the  two-dimensional  time  T,n.  That  is  r <<  r,n, 

ZD  p ZD* 

To  i t2D* 

In  one  dimension  the  detonation  moves  away  from  the 
surface  with  a velocity  VDW,  given  by 

VDW  = [2(y2-l)E0/Po]1/3  C3) 

where  y is  the  ratio  of  specific  heats  of  the  air  behind 

' the  detonation  wave,  E^  is  the  laser  flux  and  is  the 

o o 

ambient  density  in  front  of  the  wave.  The  pressure  behind 
the  wave  is 

PDW  ■ (,oVDW2/^*1)  • W 

The  surface  pressure  on  the  surface  while  the  laser  is  on 
and  the  flow  is  one-dimensional  is 

pSid  * [(Y+D/2y]2y/(Y+1)Pdw  (5) 

The  impulse  per  unit  area  at  a point  on  the  surface,  I, 
is  the  sum  of  the  contributions  to  the  impulse  transferred 
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while  the  laser  is  on  and  after  the  laser  pulse  has  terminated. 
It  is  shown  in  later  sections  of  this  report  that  a one- 
dimensional analysis  is  justified  for  the  conditions  of  the 
current  investigation  so  only  the  case  where  tq  <_  t2D  will 
be  considered. 

While  the  laser  pulse  is  on,  the  impulse  transferred  to 
the  surface  is 

t 

I = j pdt  ■ PSIDt  (0  i t < Tp)  (6) 

o 

After  the  laser  pulse  is  off  the  shock  expansion  flow  field 
relaxes  as  a planar  blast  wave  until  t2D  is  reached.  The 
pressure  is  then 

Psct)  = pSIDCTP/t)2/3  CTpit<T2D)  m 

The  impulse  delivered  to  the  surface  is  then 

T t 

1 " PS1D  + PS1D  ^ f dt/t2/3 

0 

= PS1D  Tp  t ^PS1D^PS (t) ) ^ 2 ' (0  1 t ± t2D)  (8) 


The  total  impulse  is  now  obtained  by  setting  Ps(t)  equal  to 
p , the  ambient  pressure  in  front  of  the  detonation  wave. 


' I 
< 


Then,  by  substituting  equations  3 through  5 into  equation  8 


and  by  assuming  Pgj^  >>  PQ,  it  can  be  shown  that 

i 

1/  2 

i - „■]  P°  *°TP  * {9) 

ZY  (y+ir/z  p01/z  vo 

for  x0  t2Ij,  where  VQ  is  the  ambient  speed  of  sound.  The 
impulse  delivered  to  the  surface  is,  therefore,  proportional 
to  the  speed  of  sound  in  the  ambient  gas. 

Experiments  performed  by  Pirri,  Schier  and  Northam  [Ref. 

22]  showed  the  theory  developed  above  to  predict  satis- 
factorily the  momentum  transfer  to  a surface  irradiated  with 

a high-power  pulsed  C02  laser  in  an  average  flux  regime 

2 

between  4 and  200  MW/cm  . 

Although  Pirri' s theory  presents  a comprehensive  treatment 

of  momentum  transfer  to  a surface  by  a laser  it  is  questionable 

whether  it  is  completely  applicable  to  the  conditions  used  in 

this  current  investigation.  According  to  Harrison  [Ref.  23], 

in  order  for  a laser- supported  detonation  wave  to  be  formed 

5 8 2 

laser  power  densities  of  the  order  of  10  to  10  watts/cm 
are  required  in  an  atmosphere  of  pressures  above  76  torr. 
Portions  of  the  present  investigation  were  carried  out  at 
pressures  of  the  order  of  10~5  torr.  At  these  low  pressures 
the  impulse  delivered  to  the  target  as  a result  of  shocking 
the  air  above  the  target  is  not  present.  If  this  contribution 
to  the  total  impulse  is  neglected,  the  impulse  delivered  to 


the  target  is 


pSlDt 


This  is  just  the  expression  for  impulse  per  unit  area 
delivered  to  a target  in  vacuum  reported  by  Gregg  and 
Thomas  [Ref.  24]. 

A more  suitable  equation  for  predicting  the  maximum 
pressure  developed  in  a target  is  given  by  Anderholm 
[Ref.  25].  Here  the  target  is  treated  as  a perfect  gas 
bounded  by  material  with  acoustic  impedance  Z.  The  maximum 
pressure  developed  is 


CZE„/3y1/2 


where  EQ  is  the  laser  energy  density  as  before  and  Tp  is  the 
laser  pulse  width.  The  acoustic  impedance  is  given  by 


with  p being  the  mass  density  and  VQ  the  bulk  sound  velocity 
in  the  target.  Verification  of  the  maximum  pressure 
predicted  by  this  equation  was  obtained  using  a ruby  laser 
to  generate  a pressure  of  34  kbar  in  a quartz  target.  The 
, value  for  maximum  pressure  of  29  kbar  predicted  by  equation  11 
agreed  satisfactorily  with  the  experimentally  measured  value. 


Irradiation  of  opaque  targets  can  generate  very  high 
pressures  in  the  target.  Several  experiments  have  been 
performed  by  various  authors  to  verify  this  fact.  By 
irradiating  a 0.66  mm  thick  aluminum  target  with  a pulse  of 
power  density  2000  watts/cm  from  a C02  laser,  McMordie 
and  Roberts  [Ref.  26]  measured  pulses  of  up  to  2 Bar.  As 
noted  above  pressures  in  the  kbar  range  were  measured  by 
Anderholm.  By  coating  1 mm  thick  targets  with  a water  and 
paint  coating  and  then  irradiating  the  targets  with  light 
from  a Q-switched  Neodymium: glass  laser  with  a laser 
fluence  of  25  joules/cm  , stresses  in  the  kilobar  ranges 
were  measured  by  Fox  [Ref.  27]. 

If  the  absorbing  material  is  confined  in  some  way,  as 
in  the  experiment  by  Fox  described  above,  the  mechanism 
generating  the  stresses  in  the  target  is  somewhat  changed. 

Confinement  of  the  absorbing  material  between  two  trans 
parent  films  causes  a confined  high- temperature  plasma  to 
be  produced  when  the  target  is  irradiated  by  a laser.  This 
confined  plasma  drives  a shock  wave  into  the  transparent 
media  [Ref.  28].  The  amplitude  of  this  shock  can  be  of  the 
order  of  tens  of  kilobars. 

Covering  the  absorbing  material  with  a thin  coating 
and  then  illuminating  it  with  a high-power  laser  brings  a 
different  effect  into  play.  The  thin  film  is  partially 
vaporized  as  a result  of  heating  by  thermal  conduction  from 
the  laser-heated  absorber  surface.  The  resulting  expansion 
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of  gas  and  partial  confinement  at  the  absorber/ thin  film 
interface  causes  stress  waves  to  be  generated  in  the 
absorbing  material  [Ref.  29] . Experiments  performed  by- 
Magee,  Armistead  and  Krehl  [Ref.  29]  and  by  Fox  [Ref.  27] 
have  shown  that  the  volatility  of  the  coating  has  a signifi- 
cant effect  upon  the  amplitude  of  the  stress  wave  produced. 
Magee  noted  that  as  the  curing  time  of  the  silver  paint 
used  as  a coating  for  his  targets  was  increased  the  magni- 
tude of  the  stress  wavss  produced  in  the  target  was  reduced. 
This  result  was  attributed  to  a change  in  the  transmittance 
of  the  paint  as  its  butyl  acetate  binder  evaporated.  In 
the  experiment  by  Fox  reported  above  [Ref.  26]  it  was 
observed  that  the  coating  of  the  target  with  a mere  drop  of 
water  significantly  increased  the  amplitude  of  the  stress 
wave  produced.  For  an  aluminum  target  0.17  mm  thick  it 

was  noted  that  backface  spallation  occurred  at  laser  fluences 

2 

around  310  joules/cm  when  the  target  was  water-coated. 

2 

For  an  uncoated  target  fluences  of  up  to  500  joules/cm 
produced  no  spallation.  If  the  target  was  also  coated  with 
black  paint,  the  threshold  for  backface  spallation  was 
further  reduced  to  200  joules/cm  . 


III.  EXPERIMENTAL  DESIGN 

A.  EQUIPMENT 

A Carbon  Dioxide  (COp  Transverse  Excitation  at 

Atmospheric  Pressure  (TEA)  laser  was  utilized  throughout  this 

experiment.  A description  of  the  laser  is  given  by  Strickland 

[Ref.  30],  Behrens  [Ref.  31]  and  Stevenson  [Ref.  32],  The 

' laser  output  was  reflected  off  aim  focal  length  copper 

mirror  through  a germanium  lens  onto  the  solar  cell  model. 

The  solar  cell  was,  in  turn,  mounted  in  a test  chamber  that 

could  be  maintained  either  at  atmospheric  pressure  or  at  a 
_ c 

vacuum  of  10  torr.  depending  upon  the  conditions  desired. 

The  solar  cell  model  was  mounted  on  a Sandia-type  Quartz 
Pressure  Gauge  which  enabled  the  pressure  at  the  back  face 
of  the  model  to  be  monitored  on  a Tektronix  7704  oscillo- 
scope. The  energy  output  of  the  laser  was  monitored  using 
a Ballistic  Thermopile.  A focal  spot  size  of  approximately 
1.0  cm  was  maintained  on  the  solar  cell  model  by  means  of  the 
germanium  lens.  A schematic  of  this  experimental  arrangement 
is  shown  in  Figure  3. 

1.  CC>2  TEA  Laser 

The  Naval  Postgraduate  School's  CO-  TEA  Laser  System 

L 

consists  of  a resonant  cavity,  a Tachisto  Marx-type  generator, 
a passive  electrical  component  pulse-shaping  network,  a direct 
current  power  source  and  a gas  flow  system. 
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The  resonant  cavity  consists  of  a gas-filled  lucite 


box  constructed  of  2.54  cm  thick  lucite,  101.6  cm  long, 

30.5  cm  wide  and  15.3  cm  high,  and  two  mirrors.  At  each 
end  of  the. box  a Sodium  Chloride  (NaCl)  crystal  is  mounted 
to  serve  as  a Brewster  window.  The  back  end  of  the  resonant 
cavity  is  a 10.2  cm  diameter  polished  copper  mirror  having 
a 3.0  m focal  length.  The  front  mirror  is  an  11.4  cm 
diameter  plane  parallel  germanium  mirror  which  transmits 
approximately  47%  of  the  radiation  striking  it  from  inside 
the  cavity,  thus  providing  an  outlet  for  the  laser  radiation. 

Inside  the  lucite  box  is  the  electrode  assembly  which 
consists  of  a 2.0  cm  thick  aluminum  plate  located  at  the 
bottom  of  the  box  and  a cathode  assembly  at  the  top.  The 
plate  is  cut  to  fit  the  inside  of  the  box  and  serves  as 
the  anode  of  the  system.  On  the  inside  top  of  the  box  is 
the  cathode.  It  is  composed  of  131  aluminum  blades  separated 
by  lucite  spacers  and  electrically  connected  by  two  threaded 
brass  rods.  The  trigger  assembly  consists  of  130  glass- 
encased  nichrome  wires  set  between  the  aluminum  cathode 
plates  with  each  wire  set  at  the  same  height  as  the  edge  of 
the  cathode  blade.  A gap  of  5.8  cm  is  maintained  between 
the  anode  wires  and  the  anode  plate  [Ref.  32].  Some 
characteristics  of  the  laser  are  given  in  Table  I of 
Appendix  A. 


Selective  excitation  of  the  C02  upper  laser  level 
(the  001  level)  in  a TEA  laser  is  collision  dominated.  An 


electric  discharge  across  the  gap  between  the  anode  and 
cathode  injects  energy  into  the  laser  cavity  which  creates 
a uniform  flux  of  electrons  having  an  average  energy  of  about 
one  electron  volt.  A large  number  of  inelastic  collisions 
then  take  place  between  the  CO^  molecule  in  the  ground 
state  and  the  excited  electrons.  This  causes  the  carbon 
dioxide  molecule  to  be  excited  to  various  vibrational  levels. 
Resonances  exist  in  the  carbon  dioxide  excitation  cross- 
section  for  electron  energies  of  0.3,  0.6,  and  0.9  electron 
volts.  These  resonances  coincide  with  the  vibrational 
ladder  of  energy  levels  of  the  carbon  dioxide  molecule  and 
hence  preferential  excitation  to  these  levels  occurs.  Thus 
an  electron  energy  in  the  glow  discharge  near  one  electron 
volt  excites  the  desired  energy  level  ladder.  Collisions 
between  excited  and  unexcited  carbon  dioxide  molecules  thus 
populate  the  upper  laser  level  very  efficiently.  The  high 
efficiency  arises  because  the  energy  levels  are  equally 
spaced  and  little  energy  is  lost  during  transitions  up 
or  down  the  vibrational  ladder.  This  high  excitation  effi- 
ciency, coupled  with  a quantum  efficiency  of  about  41%, 
enables  an  operational  efficiency  of  around  25%  to  be 
obtained.  (Quantum  efficiency  is  defined  as  the  ratio  of 
the  laser  photon  energy  to  the  upper  level  excitation 
energy. ) 

The  presence  of  nitrogen  ) in  the  laser  atmosphere 
provides  a complementary  excitation  process.  Nitrogen 


molecules  are  also  excited  to  a higher  vibrational  energy 
level  by  electron  impact.  The  lowest  excited  vibrational 
energy  levels  o£  nitrogen  and  the  upper  laser  level  of  the 
carbon  dioxide  molecule  are  closely  matched  so  an  efficient 
transfer  of  energy  occurs  through  collision  of  excited 
i nitrogen  molecules  with  ground  state  carbon  dioxide  molecules. 

Once  a carbon  dioxide  molecule  is  selectively 
excited  to  the  upper  laser  (001)  level,  it  can  emit  a 10.6  ym 
photon  as  it  decays,  or  is  stimulated,  to  the  lower  laser 
level  (010).  It  must  be  quickly  removed  from  this  level  to 
maintain  the  population  inversion  necessary  for  lasing.  This 
is  accomplished  by  adding  a light  gas,  such  as  Helium  (He), 
with  a collision  rate  much  larger  than  carbon  dioxide  to  the 
laser  atmosphere  to  increase  the  rate  of  depopulation  of  the 
lower  laser  level  to  the  ground  state. 

The  NPS  CC>2  laser  is  operated  with  an  atmosphere  of 
carbon  dioxide,  nitrogen  and  helium  in  a mixture  of  90/5/5 
per  cent  and  at  a flow  rate  of  about  one-half  liter  per 
second. 

Conditions  for  lasing  are  induced  by  electrical 
breakdown  between  the  cathode  and  anode  located  in  the  laser 
box.  This  produces  the  direct  electron  collision  and 
molecule-molecule  energy  transfer  process  that  results  in 
the  carbon  dioxide  population  inversion.  The  laser  is 
operated  at  atmospheric  pressure  rather  than  at  low  pressure 
in  order  to  make  use  of  the  high  molecular  density  and 
consequently  greater  energy  available. 
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The  electrical  power  required  for  excitation  of  the 
laser  is  produced  by  a three-stage  Marx  generator  supplied 
with  a positive  power  voltage  of  up  to  22  kilovolts  and  15 
milliamperes  by  an  unregulated  d.c.  power  supply.  Due  to  the 
voltage  multiplication  of  the  Marx  generator  voltages  of  up 
to  66  kilovolts  can  be  generated  across  the  laser  electrodes. 
Pulse-shaping  is  required  to  obtain  the  uniform  glow  discharge 
needed  for  gas  pumping.  This  is  accomplished  through  the  use 
of  a network  of  passive  circuit  elements  that  controls  the 
voltage  rise  time  and  decay  time  constants  of  the  input 
energy  and  the  Marx  generator.  This  network  is  discussed 
in  Ref.  31. 

2.  Solar  Cell  Model 

Reference  33  describes  the  fabrication  of  a Gallium 
Arsenide  Solar  Cell.  The  Gallium  Arsenide  (GaAs)  substrate 
is  coated  with  a continuous  platinum  (Pt)  film  approximately 
40  Angstroms  thick  to  establish  the  barrier  contact  height. 
This  platinum  film  is  then  overlaid  with  a gold  (Au)  contact 
several  thousands  of  Angstroms  thick.  Since  a gold  layer 
of  this  thickness  is  essentially  opaque  to  visible  light, 
a large  area  solar  cell  such  as  is  used  on  satellites, 
utilizes  a grid  structure  for  the  gold  contact  layer.  Such 
a grid  structure  is  illustrated  in  Figure  4.  The  entire 
solar  cell  is  then  coated  with  an  antireflective  coating  of 
silicon  dioxide  (SiC^). 

The  model  of  the  gallium  arsenide  solar  cell  con- 
structed for  the  current  investigation  was  designed  to 
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mm*** 


approximate  the  cross-section  of  the  cell  at  one  of  the 
grid  lines.  An  illustration  of  this  model  is  given  in 
Figure  1.  The  model  consisted  of  a gallium  arsenide  sub- 
strate approximately  3.0  mm  thick  and  1.0  cm  in  diameter. 

A gold  film  of  thickness  varying  up  to  5000  Angstroms  was 
vacuum  deposited  on  this  substrate.  A silicon  dioxide  film 
of  thickness  up  to  3000  Angstroms  was  then  deposited  onto 
the  gold  film.  The  model  was  then  mounted  on  a Sandia  quartz 
gauge  using  a thin  layer  of  epoxy  cement. 

Since  the  gold  layer  is  thick  enough  so  that  all 
of  the  laser  radiation  incident  upon  it  is  absorbed,  the 
momentum  coupling  between  the  gold  and  the  platinum  layers 
would  be  entirely  mechanical.  It  [Ref.  34]  has  been  stated 
that  a film  of  thickness  on  the  order  of  50  Angstroms  is 
invisible  as  far  as  mechanical  coupling  is  concerned. 
Therefore  the  thin  platinum  film  was  not  included  on  the 
solar  cell  model  used  in  the  current  investigation. 

3.  Sandia  Quartz  Gauge 

The  piezoelectric  quartz  gauges  used  in  this  experi- 
ment were  provided  by  Stanford  Research  Institute  and  were 
of  the  type  developed  by  Graham  [Ref.  5]  at  Sandia  Labora- 
tories. This  type  of  pressure  gauge  has  a resolution  on  the 
order  of  a few  nanoseconds . 

The  operation  of  the  gauge  is  based  on  the  fact 
that  when  a force  is  applied  on  a piezoelectric  crystal, 
such  as  quartz,  it  results  in  the  formation  of  a stress 
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within  the  crystal  that  distorts  the  ionic  structure  of  the 
crystal.  This  distortion  results  in  an  increase  in  the 
polarization  of  the  sample. 

It  is  well-known  from  electromagnetic  theory  that 
the  electric  displacement  vector  is  given  by 


D - eE  + P 


where  P is  the  piezoelectric  polarization,  E is  the  electric 
field  and  e is  the  dielectric  permittivity.  If  a one- 
dimensional strain  is  assumed,  the  piezoelectric  equation 
becomes 


Dx  = Td21  + eE 


where  T is  the  stress  and  d2j  is  the  piezoelectric  stress 
constant.  If  it  is  assumed  that  the  piezoelectric  polariza- 
tion is  directly  proportional  to  the  stress  with  a coefficient 
f that  is  independent  of  time  and  stress  over  a given  stress 
range,  Px  becomes 


Px  ■ foxCx) 


where  ax(x)  is  the  x-component  of  the  stress.  For  an  X-cut 


quartz  crystal  [Ref.  5] 
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d21  = 2.25*  10  coulomb/newton 


e * 4.06*  10"  farad/meter. 


For  the  static  case,  the  open  circuit  voltage  is  given  by 


[Ref.  35] 


g 3 d21/e  = °*55  &F/Tw  (volt-meter/newton) 


and  the  short-circuited  charge  for  a given  applied  force  is 


Q = d71  Fl/t 


where  F is  the  applied  force,  considered  positive  for  an 
extensional  force,  a is  the  crystal  length,  w is  the  crystal 
width  and  t is  the  crystal  thickness. 

The  development  given  above  is  adequate  for  low- 
frequency  impulses.  However,  in  the  current  investigation 
the  analysis  is  complicated  because  the  laser  pulse  produces 
a stress  wave  with  a short  wavelength. 

Graham,  et  al.  [Ref.  5]  have  analyzed  the  short- 
circuit  piezoelectric  current  induced  by  a rapidly  changing 
impulsive  load  applied  to  the  quartz.  This  load  propagates 
as  a stress  wave  along  the  x-axis  of  a quartz  disc.  The 
displacement  current  generated  by  this  stress  wave  is 


A dD/dt 
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where  A is  the  electrical  area  of  the  disc  and  D is  the 
displacement  vector. 

From  equation  14  and  assuming  one-dimensional  stress 
in  the  x-direction  it  is  seen  that 

III 
f DXW  dx  a f p(x)  dx  + e Ex(x)  dx  (19) 

0 0 0 

where  l is  the  thickness  of  the  quartz  disc. 

Assuming  a short  circuit  between  the  electrodes 
and  that  the  crystal  has  a constant  permittivity,  the  last 
integral  in  equation  (19)  becomes 

l 

J~  eEx(x)  dx  = 0 (20) 

o 

The  conductivity  of  the  quartz  is  essentially  zero, 
so 


l 

1/1  f PxCx)  d*  = Dx  * 
o 


From  equation  15  it  is  known  that 


Px  = fa  (x)  . 


(21) 
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If  it  is  assumed  that  the  quartz  is  linearly  elastic  it  can 
be  stated  that 

c(x,t)  * o (x-ugt)  (22] 

where  us  is  the  wave  propagation  velocity.  Therefore,  from 
equation  18,  we  have 

i 

ijj  = A d/dt  [1/&  j~  far(x-ust)  dx] 

o 

l 

= -fusA/£  3(ax)/3x  dx 
o 

= f us  A (o0  - ap) / Z (23] 

where  aQ  is  the  x-component  stress  at  the  input  electrode 
and  ap  is  the  x-component  at  the  rear  electrode. 

For  times  less  than  the  time  required  for  the  stress 
wave  to  travel  the  length  of  the  quartz,  ap  = 0 and 

iD  = £ A us  V*  * W 

Therefore,  for  times  less  than  the  wave  transit  time,  the 
piezoelectric  current  generated  in  the  quartz  crystal  is 
directly  proportional  to  the  stress  generated  in  the  quartz 
by  the  laser  pulse.  The  gauge  used  in  this  investigation 
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has  a time  resolution  on  the  order  of  a few  nanoseconds  and, 
using  a 50  ohm  shunt  resistor,  has  a sensitivity  of  1.09 
mV/bar.  The  gauge  is  illustrated  in  Figure  5. 

This  gauge  employed  a guard  ring  configuration  as 
described  by  Jones  [Ref.  6]  and  illustrated  in  Figure  6. 

As  is  shown,  the  inner  portion  of  the  X-cut  quartz  disc  is 
electrically  isolated  from  the  outer  portion  by  separating 
the  vapor-coated  electrode  into  two  regions.  This  guard 
ring  configuration  eliminates  field  distortions.  At  the  outer 
edge  of  the  disc,  the  discontinuity  in  electric  potential 
and  dielectric  permitivity  of  the  quartz  produces  electric 
field  fringing  analogous  to  that  found  in  parallel-plate 
capacitors.  The  guard  ring  restricts  observation  to  the 
central  region  of  the  quartz  disc  where  there  is  negligible 
fringing. 

This  guard  ring  also  aids  in  the  elimination  of 
distortions  caused  by  unloading  waves.  As  the  stress  wave 
propagates  through  the  disc  in  the  axial  direction,  boundary 
conditions  generate  shear  and  dilational  waves  immediately 
behind  the  wavefront  at  the  lateral  edge  of  the  disc.  These 
unloading  waves  then  propagate  laterally  inward  from  the 
edge  and  the  crystal  is  no  longer  in  a state  of  one -dimensional 
strain.  The  central  region  of  the  disc  will  be  one-dimensional 
only  for  the  first  wave  transit  time  if  the  width  of  the 
center  electrode  is  such  that  the  unloading  wave  does  not 
reach  the  central  region  during  the  first  wave  transit  time. 


For  proper  monitoring  of  the  pressure  pulse,  the 


propagation  time  of  the  pressure  pulse  in  the  quartz  must 
be  larger  than  the  half-width  of  the  laser  pulse.  For  the 
250  ns  pulse  width  of  the  C02  TEA  laser  and  a speed  of  sound 
in  quartz  of  5700  meters  per  second,  this  requires  that  the 
gauges  used  in  the  CC>2  laser  must  have  quartz  discs  at  least 
1.44  mm  thick.  The  gauges  used  had  quartz  discs  3.0  mm 
thick  so  this  requirement  was  satisfied. 

Care  also  had  to  be  taken  that  a pulse  reflected  at 
the  rear  of  the  gallium  arsenide  target  sample  did  not  reach 
the  gauge  during  the  loading  phase.  This  required  that  the 
target  sample  thickness,  dg,  had  to  satisfy  [Ref.  29] 

ds  2.  c,  V2  (2sj 

where  t^  is  the  propagation  time  of  a pressure  (sound)  pulse  in 
the  quartz  plate  and  vg  is  the  speed  of  sound  in  the  gallium 
arsenide.  For  a sound  speed  of  5000  meters  per  second  in 
gallium  arsenide  and  the  500  ns  transit  time  in  the  quartz 
disc  calculated  from  the  data  for  quartz  given  above,  this 
/required  that  dg  be  greater  than  1.25  mm.  The  3.0  mm  thickness 
of  the  gallium  arsenide  samples  used  satisfied  this  requirement. 

A final  requirement  for  proper  operation  of  the  quartz 
gauge  was  that  the  stress  be  essentially  uniaxial.  To  ensure 
this  a laser  spot  diameter  of  slightly  greater  than  the  1.0 
cm  diameter  of  the  gallium  arsenide  sample  was  used 
throughout  the  experiment. 


i 
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. Ballistic  Thermopile 

The  total  energy  output  of  the  C02  TEA  laser  was 
monitored  by  directing  the  laser  beam  into  a ballistic 
thermopile.  The  ballistic  thermopile  is  a calorimetric 
device  that  measures  the  temperature  rise  due  to  absorbed 
radiation.  It  is  composed  of  two  nickel-plated  cones 
(receiver  and  reference)  and  two  series-connected  iron- 
constantan  thermocouples.  These  thermocouples  are  situated 
so  that  the  hot  junctions  are  attached  to  the  receiver  cone 
and  the  cold  junctions  to  the  reference  cone.  Radiant 
energy  is  directed  into  the  receiver  cone  where  it  is  almost 
totally  absorbed.  The  temperature  difference  between  the 
* receiver  cone  and  reference  cone  generates  an  emf , which  is 
monitored  by  an  externally  connected  microvoltmeter.  The 
peak  meter  reading  is  linearly  proportional  to  the  total  input 
energy  with  a constant  of  proportionality  equal  to  44.9  pV/ 
joule.  Due  to  the  high  reflectivity  in  the  cone,  the  pol- 
ished surface  resists  the  destructive  effects  of  the  high 
peak  power  pulses  found  in  laser  plasma  research.  Some 
characteristics  of  the  HADRON  Ballistic  Thermopile  used  in 
this  investigation  are  given  in  Table  II  of  Appendix  A. 

5.  Test  Chamber  for  C02  Laser  System 

The  test  chamber  utilized  was  constructed  of  a 
rectangular  aluminum  block  11.4  cm  long,  8.9  cm  wide  and 
10.0  cm  high  through  which  a 6.4  cm  diameter  hole  had  been 
drilled  from  the  top  to  the  bottom.  Five  access  ports  were 
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provided,  a 4.0  cm  diameter  port  in  each  side  and  a 6.4  cm 
diameter  port  on  the  top  of  the  chamber.  An  overhead  view 
of  the  test  chamber  is  shown  in  Figure  7. 

Laser  light  was  introduced  into  the  chamber  through 
Port  A.  Ports  B and  C were  covered  with  plexiglass  windows 
to  enable  the  system  to  be  placed  under  vacuum  by  means  of 
a vacuum  system  connected  at  Port  D.  The  test  chamber 
was  mounted  on  an  aluminum  bed  plate  and  that  was  adjustable 
in  height  by  means  of  bolts  mounted  at  each  corner  of  the 
plate . 

B . THEORY 

1.  Determination  of  Peak  Pressure  in  the  Gallium 

Arsenide  Substrate 

Equation  10  provided  a method  whereby  the  maximum 
pressure  developed  in  a transparent  material  by  laser 
irradiation  can  be  calculated.  However,  the  problem 
presented  by  the  generation  of  a pressure  pulse  in  the  GaAs 
layer  of  the  solar  cell  model  used  in  the  current  investi- 
gation is  complicated  by  the  fact  that,  although  GaAs  is 
essentially  transparent  to  the  10.6  ym  radiation  used,  the 
gold  and  silicon  dioxide  coatings  are  strongly  absorbing 
at  this  wavelength.  Also,  the  thickness  of  these  two  films 
was  such  that  the  incident  laser  radiation  was  totally 
absorbed  prior  to  reaching  the  GaAs  substrate.  Hence  the 
problem  of  predicting  the  magnitude  of  the  pressure  pulse 
in  the  GaAs  was  reduced  to  that  of  predicting  the  pressure 
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in  the  gold-silicon  dioxide  combination  layer  and  then 
computing  a correction  factor  for  the  mechanical  coupling 
of  the  pressure  pulse  across  the  gold-gallium  arsenide 
interface. 

By  consideration  of  the  conservation  of  momentum 
across  the  interface  between  two  materials  of  different 
acoustical  properties  Krehl,  et  al.  [Ref.  19]  developed  an 
equation  for  this  correction  factor.  Employing  this  result, 
the  magnitude  of  the  pressure  in  the  gallium  arsenide 
transmitted  by  the  gold  layer  is 

PAu  " PGaAs  ^vGaAspGaAs  * vAu°Au)/2  vGaAspGaAs  <26> 

where  p is  the  initial  density  and  v is  the  elastic  wave 
velocity  in  the  medium.  Using  the  values  of  p and  v for 
gallium  arsenide  and  gold  given  in  Table  III  of  Appendix  A 
we  obtain 


PGaAs/PAu  ' °-43  • (2?) 

A similar  correction  must  be  made  to  the  pressure  measured 
by  the  quartz  gauge  to  obtain  the  pressure  in  the  gallium 
? 'senide.  In  this  case  the  correction  factor  was  found  to 
be 

PGaAs/PSi02  = 1.57  (27a) 
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Using  equations  27  and  27a  it  is  possible  to  obtain  the  true 
pressure  generated  in  the  gallium  arsenide  if  the  pressure 
in  the  gold  layer  or  the  pressure  measured  by  the  quartz 
gauge  is  known. 

The  problem  of  predicting  the  pressure  generated  in 
the  gallium  arsenide  by  the  laser  pulse  is  now  reduced  to 
i predicting  the  pressure  in  the  gold  layer.  From  equation  11 
it  is  known  that 

PAu  ■ (2AuVV1/2  • C28) 

But  since  the  silicon  dioxide  layer  covering  the  gold  layer 
is  an  absorber  of  10.6  um  radiation  not  all  of  the  laser 
energy  incident  upon  the  silicon  dioxide  should  be  trans- 
mitted to  the  gold  layer.  However,  from  Table  III  of 
Appendix  A,  it  is  seen  that  the  absorption  coefficient  of 
silicon  dioxide  at  10.6  pm  is  3.0’103  cm  1.  For  the  silicon 
dioxide  films  with  thicknesses  on  the  order  of  thousands  of 
Angstroms,  such  as  those  used  in  the  current  investigation, 
the  energy  fraction  transmitted  through  the  silicon  dioxide 
i layer  is  about  0.98.  So  virtually  all  of  the  incident 
energy  is  transmitted  to  the  gold  layer  by  the  silicon  dioxide 
layer  and  it  can  be  stated  that 


E 


Au 
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(29) 
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where  A is  the  absorptance  of  the  Au-Si02  film  system. 
Equation  28  then  becomes 

PAU  ' <ZAu  A • (30) 

So  if  the  laser  energy  density  Eq  and  the  absorptance  A 
are  known  the  pressure  in  the  gallium  arsenide  can  be 
calculated  using  equations  27  and  30. 

2 . Derivations  of  the  Absorptance 

Since  both  the  gold  and  the  silicon  dioxide  layers 
are  absorbers  of  10.6  ym  radiation,  the  problem  of  finding 
the  absorptance  of  the  system  is  that  of  finding  the 
absorptance  of  a single  absorbing  layer  (the  silicon  dioxide) 
on  an  absorbing  substrate  (the  gold).  This  problem  has  been 
solved  by  Heavens  [Ref.  36]. 

The  energy  that  is  coupled  into  the  system  consists 
of  the  laser  energy  that  is  absorbed  by  the  silicon  dioxide 
layer  and  the  laser  energy  that  is  transmitted  through  the 
silicon  dioxide  layer  and  is  then  absorbed  by  the  gold 
layer.  If  the  absorptance  is  defined  as  the  ratio  of  the 
energy  coupled  into  the  system  by  both  of  these  means  to 
the  incident  energy,  and  the  reflectance  of  the  system  is 
defined  as  the  ratio  of  the  energy  reflected  to  the  energy 
incident,  then 


A = 1 - R 

where  R is  the  reflectance. 


(31) 
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From  Heavens,  the  reflectance  of  a system  of  thin 


films  such  as  illustrated  in  Figure  6 is 

R = cc*/aa*  (32) 

where  c and  a are  obtained  from  the  matrix  product  of  the 
matrices  of  Fresnel  coefficients  for  the  system.  That  is, 


(C1)(C2)....(Cn+1)  = /a 


(33) 


In  general,  for  an  absorbing  medium,  the  elements  of  the 
matrix  of  Fresnel  coefficients  are  all  complex  and  are 
given  by 
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in  which  nR  and  n^  are  the  real  and  imaginary  parts  of  the 
complex  index  of  refraction  for  the  material  composing  the 
layer  and  d is  the  thickness  of  the  layer. 

For  normal  incidence  the  quantities  of  equation  34 
are  defined  by 


U.  - 

m-1 

pm  = e cos  y , 

rm  m-l 


q = e sin  y , 
^m  m-1 


r = e m'  (g  cos  y„  ■>  * h sin  Ym  ) 
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A double  suffix  notation  can  be  adopted  to  denote 
the  elements  of  the  product  matrices.  Thus  the  elements  of 
the  single  layer  system  under  consideration  are 


(cx)(c2)  = 


fp12  + iql2  ri2  + 1S12 


^12  + iu12  V12  + iw 
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where 
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u12  h2  + hlP2  " glq2  (50) 

and  the  other  terms  of  the  product  matrix  are  obtained  in 
the  same  fashion  as  equations  47  through  50  by  matrix 
multiplication  of  the  individual  matrices  (c  ) . From 
equation  32  and  from  the  quantities  defined  in  equations 
35  through  50,  the  reflectance  of  the  system  is 

R - <‘12  * u12)/tpl2  * «12J  (51> 

The  absorptance  of  the  system  can  then  be  calculated  from 
equation  31. 

In  general  the  quantities  n^  and  n^  needed  to  calcu- 
late the  terms  of  equation  51  are  wavelength  dependent. 

Values  for  these  quantities  are  tabulated  in  Table  III  of 
Appendix  A ^or  the  wavelength  of  10.6  qm  used  in  the 
current  investigation. 
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IV.  EXPERIMENTAL  PROCEDURE 

A.  MEASUREMENT  OF  PRESSURE  PULSES 

The  magnitude  of  the  pressure  pulse  at  the  backface 
of  the  gallium  arsenide  was  monitored  using  the  Sandia 
Quartz  Gauge  previously  described.  The  output  of  the 
gauge  was  displayed  on  the  screen  of  a Tektronix  7704 
Oscilloscope  and  recorded  photographically.  Triggering  for 
this  oscilloscope  was  obtained  from  a Tektronix  T-555 
oscilloscope  utilizing  its  delayed  trigger  capability. 

The  T-555  oscilloscope  was,  in  turn,  triggered  by  the 
current  output  of  the  Marx  bank  of  the  Carbon  Dioxide  laser. 
Delays  of  up  to  1.0  microsecond  were  used  to  allow  for  the 
delay  between  firing  the  Marx  bank  and  initiation  of  the 
laser  pulse.  A schematic  of  the  circuitry  used  is  shown 
in  Figure  8. 


B.  MEASUREMENT  OF  ENERGY  DELIVERED  TO  TARGET 

Prior  to  initiation  of  the  experimental  data  runs  the 
Carbon  Dioxide  laser  was  aligned  so  that  its  energy  output 
at  a given  applied  voltage  varied  less  than  ten  per  cent. 
After  the  data  run  was  completed  the  energy  output  was 
measured  using  the  ballistic  thermopile  to  verify  that  no 
significant  change  in  the  laser  output  occurred.  The  experi- 
mental apparatus  was  such  that  continuous  monitoring  of  the 
laser  output  was  not  practical. 
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Since  the  laser  beam  was  required  to  pass  through  various 
optical  components  prior  to  irradiating  the  solar  cell  model 
not  all  of  the  energy  of  the  laser  pulse  reached  the  target; 
some  was  absorbed  by  the  various  optical  components. 

The  principal  components  that  absorbed  the  laser  energy 
were  the  germanium  lens  used  to  focus  the  laser  beam  and, 
in  the  case  of  the  vacuum  data  runs,  the  window  in  the  test 
chamber  through  which  the  laser  light  was  admitted  to  the 
chamber. 

An  experiment  was  performed  in  which  the  energy  output 
of  the  laser  was  measured  with  the  ballistic  thermopile 
both  before  and  after  the  laser  beam  passed  through  the 
germanium  lens.  The  lens  was  found  to  transmit  about  90% 
of  the  laser  energy. 

A similar  experiment  was  performed  on  the  germanium 
window  used  as  an  entrance  window  for  the  vacuum  data  runs. 
However,  as  is  described  in  the  Experimental  Results  section 
of  this  report,  the  energy  transmitted  through  the  germanium 
window  was  a function  of  the  number  of  laser  shots.  A trans- 


mittance curve  was  obtained  for  this  window  and  is  shown 
' in  Figure  9.  A correction  factor  was  then  applied  to  the 
laser  energy  for  the  value  appropriate  to  the  laser  shot. 

From  the  above  experiments  a correction  factor  of  0.9 
was  applied  to  the  laser  energy  to  account  for  the  absorp- 
tance  of  the  lens  and  a correction  factor  corresponding 
to  the  transmittance  of  the  window  for  a given  laser  shot 
was  applied  to  account  for  the  absorptance  of  the  window. 
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. EXPERIMENTAL  RESULTS  AND  ANALYSIS 


I 


A.  ENTRANCE  WINDOW  DAMAGE 

An  unexpected  result  that  arose  in  the  performance  of 

the  experiment  was  the  inability  to  find  a material  for  the 

entrance  window  to  the  te*t  chamber  that  was  capable  of 

7 2 

withstanding  the  power  density  of  10  watts/cm  needed  to 
obtain  the  1.0  cm  spot  diameter  required  on  the  solar  cell 
model.  Lack  of  a window  that  could  withstand  more  than  a 
few  shots  at  this  power  density  severely  limited  the  amount 
and  quality  of  the  data  obtained  in  the  vacuum  data  runs. 

Four  different  materials  were  tried  in  an  attempt  to  find 
a suitable  material:  Sodium  Chloride,  Potassium  Bromide, 

Zinc  Selenide,  both  with  and  without  an  anti-reflective 
coating,  and  anti-reflective  coated  Germanium. 

Of  these  materials  the  Potassium  Bromide  window  showed 
damage  in  the  form  of  a burn  spot  on  the  first  shot  and  the 
Sodium  Chloride  and  the  anti-reflective  coated  Zinc  Selenide 
windows  showed  similar  damage  after  about  three  shots. 

Closer  examination  of  the  Zinc  Selenide  window  showed 
indications  of  possible  internal  damage  and  spallation  on 
its  back  face. 

Since  the  anti-reflective  coating  on  the  window  was 
damaged,  it  was  decided  to  remove  it  and  determine  if  the 
uncoated  surface  was  better  able  to  pass  the  laser  beam 
without  damage.  It  was  found  that  this  was  not  the  case. 
After  three  shots,  it  was  observed  that  further  surface 
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damage  was  occurring  and  experimentation  with  this  window 
was  terminated. 


Finally,  an  anti-reflective  coated  germanium  window  was 
obtained  and  was  found  to  survive  five  shots  without  any 
apparent  damage.  A vacuum  data  run  was  then  initiated  using 
this  window.  However,  after  a few  more  shots  surface  damage 
was  visible  on  the  window  and  the  run  terminated. 

A subsequent  data  run  was  made  by  irradiating  a dif- 
ferent portion  of  the  window  surface.  In  this  run  surface 
damage  was  observed  after  five  shots. 

Pictures  of  the  damage  caused  to  the  Sodium  Chloride 
and  Germanium  windows  are  included  as  Figures  10  and  11. 

The  primary  cause  of  the  damage  observed  in  the  window 
materials  examined  is  felt  to  be  that  the  incident  power 
density  exceeded  the  intrinsic  threshold  for  damage  of  the 
material  coupled  with  damage  due  to  possible  inclusions  in 
the  materials 

The  possibility  that  inclusions  in  the  surface  layer  of 
the  Sodium  Chloride  and  Potassium  Bromide  windows  were 
responsible  for  the  damage  to  these  windows  is  extremely 
> plausible.  These  windows  were  polished  by  the  experimenter 
and  it  is  quite  likely  that  small  pieces  of  the  polishing 
compound  could  have  been  imbedded  in  the  surface  of  the 
windows.  These  particles  could  then  have  acted  as  strong 
absorbers  for  the  incident  10.6  ym  radiation. 

The  most  probable  mechanism  to  explain  the  damage 
observed  in  the  Zinc  Selenide  and  Germanium  windows  is 
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exceeding  the  intrinsic  damage  threshold  for  the  anti- 
reflective  coatings  of  these  windows.  While  no  information 
for  the  coating  of  the  Zinc  Selenide  window  was  available, 
information  given  by  the  manufacturer  of  the  Germanium  window 
indicated  that  the  damage  threshold  for  the  coating  on  this 
window  was  about  100  megawatts/cm  . Measurement  of  the  burn 
spot  on  the  Germanium  window  indicated  that  it  had  a diameter 
of  about  0.5  cm;  this  corresponds  to  a power  density  of 
about  7000  megawatts/cm  , well  in  excess  of  the  damage 
threshold  of  the  coating. 

B.  PRESSURE  PULSE  MEASUREMENTS  IN  VACUUM 

Two  data  runs  were  made  under  a vacuum  of  approximately 
» 10  ^ torr.  using  the  germanium  entrance  window,  a laser 

energy  of  5.25  joules,  and  two  different  sets  of  thicknesses 
for  the  silicon  dioxide  and  gold  coatings  of  the  solar  cell 
model.  Tables  IV  and  V of  Appendix  A list  the  results  for 
each  set  of  film  thicknesses.  A graph  of  the  natural 
logarithm  of  the  peak  pressure  measured  in  the  gallium 
arsenide  versus  shot  number  of  the  laser  for  each  set  of 
coating  thicknesses  is  displayed  in  Figures  12  and  14. 

As  can  be  seen  on  these  graphs,  the  peak  pressure  in 
each  case  steadily  decreased  from  a high  initial  value 
to  about  100  bars  as  the  shot  number  increased.  The  dura- 
tion of  the  pressure  pulse  in  both  cases  remained  approxi- 
mately constant  at  about  400  ns  until  the  last  two  shots 
when  it  shortened  to  around  250  ns.  Figures  16  and  17 
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display  early  and  final  pressure  pulses  for  the  system  of 
a 1000  Angstrom  layer  of  silicon  dioxide  on  a 1000  Angstrom 
layer  of  gold  (1000/1000  system).  Similar  pulses  for  a 
1000/2000  system  are  shown  in  Figures  18  and  19. 

Observations  of  each  gauge  after  completion  of  the  data 
run  revealed  that  in  both  cases  the  layers  of  silicon 
dioxide  and  gold  had  been  removed  from  most  of  the  surface 
of  the  gallium  arsenide  substrate.  Since  the  duration  of  the 
data  run  for  the  thicker  of  the  two  systems  (the  1000/2000 
system)  was  the  shorter  (10  shots)  it  is  concluded  that  it 
required  an  equal  number  of  shots  to  remove  the  coatings 
from  the  1000/1000  system. 

C.  PRESSURE  PULSE  MEASUREMENTS  AT  ATMOSPHERIC  PRESSURE 

Three  data  runs  were  made  at  atmospheric  pressure  using 
the  following  sets  of  gold  and  silicon  dioxide  layer  thick- 
nesses : 

1.  1000/2000 

2.  2000/3000 

3.  0/5000 

In  the  above  sets  the  thickness  of  the  silicon  dioxide  layer 
in  Angstroms  is  followed  by  the  thickness  of  the  gold 
layer  in  Angstroms. 

In  all  cases  a laser  energy  output  prior  to  any  optical 
components  of  5.25  joules  was  used  and  a spot  diameter  of 
1.3  cm  was  maintained  on  the  solar  cell  model.  The  data 
from  these  runs  are  tabulated  in  Tables  VI  and  VII  of 
Appendix  A. 
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In  the  1000/2000  system  the  peak  pressure  pulse  in  the 

s 

gallium  arsenide  was  initially  14.0  bars  and  decreased  to 
a value  of  9.0  bars  after  13  shots  and  remained  nearly 
constant  at  this  value  for  the  remaining  six  shots  of  the 
run.  The  duration  of  the  pressure  pulses  measured  remained 
nearly  constant  at  1.0  usee.  After  six  shots,  damage  to 
the  coatings  was  observed  in  the  form  of  total  removal  of 
the  silicon  dioxide  and  gold  layers  at  various  points  on 
the  surface  of  the  solar  cell  model.  Figure  20  displays 
a typical  pulse  for  this  system.  Figure  21  displays  a 
graph  of  the  natural  logarithm  of  peak  pressure  versus  shot 
number  for  this  data. 

The  second  system  examined  at  atmospheric  pressure 
consisted  of  thicker  coatings.  In  this  case  a 2000  Angstrom 
layer  of  silicon  dioxide  was  deposited  onto  a 3000  Angstrom 
layer  of  gold.  The  results  of  the  fourteen  shots  taken  at 
this  2000/3000  system  did  not  display  the  decrease  in  peak 
pressure  observed  in  the  previously  described  results.  For 
this  system  the  peak  pressure  measured  remained  nearly 
constant  with  an  average  value  of  8.0  bars  and  displayed  a 
mean  pulse  duration  of  approximately  600  nsec. 

The  greater  thickness  of  the  silicon  dioxide  coating 
allowed  an  observation  to  be  made  to  obtain  an  estimate  of 
the  number  of  shots  required  to  at  least  partially  remove 
- the  silicon  dioxide  layer.  It  is  estimated  that  the  silicon 
dioxide  layer  was  almost  totally  removed  from  the  central 
portion  of  the  surface  of  the  solar  cell  model  after  only 
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one  shot.  After  a total  of  six  shots,  gallium  arsenide 
was  visible  at  the  center  of  the  surface  indicating  that 
both  the  silicon  dioxide  and  gold  layers  were  removed.  A 
typical  pressure  pulse  measured  for  this  system  is  shown  in 
Figure  22.  Figure  24  shows  the  surface  of  the  solar  cell 
model  upon  completion  of  fourteen  shots. 

A final  series  of  shots  was  run  on  a target  upon  which 
a 5000  Angstrom  layer  of  gold  was  deposited.  No  layer  of 
silicon  dioxide  was  deposited  on  this  solar  cell  model  in 
order  that  the  effects  of  a thick  layer  of  gold  upon  the 
pressure  pulse  generated  could  be  examined.  An  initial  peak 
pressure  of  approximately  six  bars  was  measured  with  a 
duration  of  approximately  2.0  ysec.  This  magnitude  of  the 
pressure  pulse  persisted  for  four  shots  and  then  the  pulse 
was  lost  in  the  noise  displayed  on  the  screen  of  the 
oscilloscope.  The  run  was  terminated  at  eight  shots  with 
no  damage  to  the  gold  surface  apparent.  Figure  23  shows  a 
typical  pulse  from  this  series. 

D.  COMPUTER  MODELING  OF  PEAK  PRESSURE  PULSE 

A computer  code  was  developed  to  predict  both  the  peak 
pressure  in  vacuum  and  in  atmosphere  for  a given  set  of 
silicon  dioxide  and  gold  layer  thicknesses.  For  irradiation 
in  vacuum  the  plasma  tamping  model  presented  by  equation  30 
was  used  and  for  irradiation  at  atmospheric  pressure  the 
laser-supported  detonation  wave  model  developed  by  Pirri  and 
expressed  by  equation  10  was  used. 
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In  both  models  the  resultant  pressure  is  dependent  upon 
the  energy  density  delivered  to  the  target.  From  equations 

29  and  31  it  can  be  seen  that  this  energy  density  is  at  least 
partially  dependent  upon  the  reflectance  of  the  target 
irradiated.  Consequently,  the  computer  code  was  designed 

to  first  calculate  the  reflectance  of  the  solar  cell  model 
from  equation  32,  and  from  this  to  obtain  the  absorptance 
of  the  system  from  equation  31.  The  actual  energy  density 
delivered  to  the  target  was  then  calculated  from  equation  29. 
This  resultant  energy  was  then  applied  to  equations  10  and 

30  to  obtain  the  expected  magnitudes  of  the  pressure  pulses 
in  vacuum  and  at  atmospheric  pressure  respectively.  For 
the  case  of  the  atmospheric  measurements  the  duration  of 
the  pulse  was  measured  from  the  photographic  records  of 
the  data  and  this  value  was  then  assigned  to  the  variable 

t in  equation  10. 

The  code  was  also  designed  to  examine  the  effects  of 
reducing  the  thickness  of  the  silicon  dioxide  and  gold 
layers  upon  the  peak  pressure  expected.  f'\nce  the  obvious 
dependence  of  the  predicted  pressure  on  the  thickness  of 
the  coating  layers  is  through  the  exponential  terms  of  the 
absorptance,  this  was  accomplished  by  stripping  off  the 
silicon  dioxide  and  gold  films  in  50-Angstrom  increments 
and  monitoring  the  effect  this  had  on  the  peak  pressures 
predicted. 

The  results  of  running  this  computer  code  for  the 
thicknesses  of  silicon  dioxide  and  gold  used  for  the  data 
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runs  of  the  current  investigation  are  summarized  in  Table  VIII 
of  Appendix  A.  The  computer  code  itself  is  given  in  Appendix  F. 

As  is  shown  in  this  table,  in  all  cases  where  a silicon 
dioxide  layer  was  present  the  initial  reflectance  of  the 
system  was  almost  99%.  At  the  point  where  the  entire 
silicon  dioxide  layer  is  removed,  the  reflectance  is  only 
reduced  by  about  1%.  The  value  of  the  reflectance  at  this 
point  corresponds  to  that  of  the  uncoated  gold  layer,  that 
is  about  9%. 

E.  ANALYSIS  OF  EXPERIMENTAL  RESULTS 

1 . Pressure  Measurements  in  Vacuum 

The  data  depicted  in  Figures  12  and  14  was  fitted 
to  a straight  line  using  a linear  regression  computer  code. 

This  was  done  to  test  the  hypothesis  that  the  peak  pressure 
measured  had  a functional  dependence  upon  shot  number  of 
the  form 

P = PQ  ekx.  (52) 

where  k is  some  constant  dependent  upon  the  initial  thicknesses 
of  the  components  of  the  film  system  and  x is  the  shot  number. 
The  linear  regression  analysis  yielded  a value  for  k,  the 
peak  pressure  to  be  expected  after  the  first  shot  of  the 
laser  and  a correlation  factor  that  was  a measure  of  the 
degree  to  which  the  data  fitted  the  assumed  straight  line. 

Since  the  data  obtained  in  the  vacuum  runs  was 
complicated  by  the  fact  that  the  entrance  window  of  the 


chamber  absorbed  some  portion  of  the  laser  energy  that 
varied  from  shot  to  shot,  the  measured  pressure  was 
corrected  for  the  fraction  of  energy  absorbed  by  the  window. 
This  was  accomplished  by  dividing  the  pressure  for  a given 

shot  number  by  the  square  root  of  the  transmittance  of  the 

1/2 

germanium  window  for  that  shot  number.  The  factor  T ' 

was  chosen  since  the  energy  transmitted  through  the  window 

is  TEq  where  E^  is  the  energy  incident  on  the  window.  The 

plasma  tamping  model  used  for  predicting  the  pressure  in 

vacuum  is  expressed  by  equation  30  wherein  the  energy 

density  enters  to  the  one-half  power;  hence  the  correction 
1/2 

factor  T is  applied. 

Dividing  by  this  factor  has  the  effect  of  normalizing 
all  the  pressures  to  the  values  that  would  have  been  measured 
had  the  window  not  been  present.  A check  to  see  if  the  only 
effect  the  absorption  of  energy  by  the  window  was  to  lower 
the  pressure  measured  by  comparing  the  value  of  k calcu- 
lated by  fitting  a straight  line  to  this  "corrected”  data 
with  that  obtained  by  analysis  of  the  uncorrected  data. 
Figures  12  and  14  display  a graph  of  these  corrected 
pressures  versus  shot  number.  If  the  values  for  k were 
equal  or  nearly  so  then  the  dependence  of  the  peak  pres- 
sure upon  film  thickness  was  independent  of  the  energy 
absorption  by  the  window. 

For  the  1000/1000  system  the  analysis  yielded  a value 
of  k equal  to  -0.14  before  the  correction  for  the  window 
was  applied;  after  the  correction  was  applied  k was  equal 
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to  -0.13.  For  the  1000/2000  system  the  uncorrected  value  of 
k was  -0.15  and  after  correction  was  -0.15.  Based  on  the 
good  agreement  between  the  values  of  k before  and  after 
the  correction  was  applied,  the  assumption  that  the  energy 
- absorption  by  the  window  was  not  dependent  upon  the  thickness 
of  the  films  was  considered  justified. 

To  determine  the  validity  of  the  hypothesis  of  the 
functional  relationship  between  shot  number  and  peak  pressure 
indicated  by  equation  52,  the  correlation  factor  calculated 
by  the  linear  regression  code  was  examined.  In  the  case 
of  the  1000/1000  system  a correlation  of  0.79  was  calculated. 
For  the  1000/2000  system  the  correlation  factor  was  0.89. 
Since  a correlation  of  1.0  indicates  a perfect  straight 
line,  the  values  obtained  indicated  that  the  functional 
dependence  assumed  in  equation  52  was  possibly  justified. 

However,  the  correlation  value  of  0.79  calculated 
for  the  1000/1000  system  seemed  to  indicate  that  the 
functional  relationship  between  peak  pressure  and  shot 
number  of  equation  52  was  not  entirely  correct.  Closer 
examination  of  the  data  for  the  1000/1000  system  displayed 
in  Figure  11  indicated  that  the  peak  pressure  measured 
remained  relatively  constant  from  shot  10  on.  Consequently, 
a functional  relationship  of  the  form 

— k ,x 

P = P + PQeK  (53) 
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' where  k'  was  assumed  to  be  thickness  dependent  was  assumed 
for  the  first  nine  shots  of  the  1000/1000  and  1000/2000 
systems.  After  nine  shots  it  was  presupposed  that  the 

I 

-pressure  saturated  to  F,  the  average  pressure  measured  for 
shots  10  through  16  of  the  100/100  system.  Since  the  data 
for  the  1000/2000  system  covered  only  ten  shots  it  was  not 
possible  to  determine  a corresponding  value  of  P for  that 
system  and  the  value  for  the  1000/1000  system  was  also  used 
for  the  1000/2000  system. 

To  verify  the  relationship  suggested  by  equation  53 
the  mean  value  of  the  pressures  measured  for  shots  ten  through 
sixteen  of  the  1000/1000  system  was  calculated  and  found  to 
be  215  bar.  This  value  was  then  subtracted  from  P,  the 
value  of  peak  pressure  after  correction  for  the  absorptance 
of  the  germanium  window,  corresponding  to  each  shot  number 
for  the  results  of  the  vacuum  data  runs.  The  results  of 
this  subtraction  are  displayed  in  Tables  IX  and  X of  Appendix 
A.  A graph  of  the  natural  logarithm  of  the  resultant  pressure 
P'  versus  shot  number  was  then  made.  These  graphs  are  shown 
as  Figures  13  and  15.  A linear  regression  analysis  was 
performed  on  these  data  to  test  the  validity  of  the 
relationship 

f 

P»  « (P  - P)  - Poek'X  (54) 

suggested  by  equation  53  for  the  first  nine  shots. 


For  the  1000/1000  system  a value  of  k'  of  -0.44 


and  a correlation  of  0.91  were  obtained.  The  corresponding 
results  for  the  1000/2000  system  were  k'  = -0.22  and  a 
correlation  of  0.85.  Comparing  these  correlations  to  those 
reported  above  for  the  verification  of  equation  52  it  is 
seen  that,  in  the  case  of  the  1000/1000  system,  the  correla- 
tion increased  significantly.  For  the  other  system  considered, 
the  1000/2000  system,  the  correlation  did  not  change 
appreciably. 

Viewing  the  fact  that  equation  54  yielded  a higher 
correlation  for  the  1000/1000  system  than  did  equation  52 
and  that  the  results  for  the  1000/2000  system  did  not 
appreciably  change,  it  was  concluded  that  the  relationship 
between  shot  number  and  pressure  generated  is  of  the  form 
of  equation  54.  That  is,  a "saturated"  exponential  for  the 
first  nine  shots  and  a constant  value  thereafter.  The 
saturation  is  felt  to  occur  when  the  gold  and  silicon 
dioxide  coating  are  entirely  removed. 

An  assumption  was  made  that  the  value  of  k in 
equation  52  and  k'  in  equation  54  are  dependent  upon  the 
initial  thicknesses  of  the  silicon  dioxide  and  gold  films 
of  the  solar  cell  model.  An  attempt  to  at  least  partially 
justify  this  assumption  was  made  using  the  previously 
described  computer  code. 

If  the  relationship  between  peak  pressure  and  film 
thickness  exists  and  is  of  the  form  of  equation  54,  that  is, 
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a saturated  exponential  dependent  upon  the  initial  thicknesses 
of  the  gold  and  silicon  dioxide  films,  the  dependence  will 
appear  in  the  absorptance  of  the  system  since  that  is  the 
only  place  where  the  film  thicknesses  appear.  If  there  is 
a relationship  between  the  constant  k'  in  equation  54  and 
the  initial  thicknesses  of  the  films  then  the  pressure 
predicted  by  the  plasma  tamping  model  should  equal  that 
measured  after  the  first  shot  of  the  laser. 

The  initial  pressures  measured  for  the  film  systems 
used  (as  predicted  after  the  data  had  been  analyzed  by  the 
linear  regression)  were: 

1.  1000/1000  system:  1706  bar. 

2.  1000/2000  system:  753  bar. 

The  values  predicted  using  equation  30  were: 

1.  1000/1000  system:  556  bar. 

2.  1000/2000  system:  556  bar. 

The  predicted  value  for  each  case  is  respectively  33%  and 
74%  of  the  measured  value.  Hence,  it  is  concluded  that 
there  is  at  least  some  justification  to  the  assumed  rela- 
tionship between  film  thickness  and  the  proportionality 
constant  k'  of  equation  54.  The  existence  of  this  rela- 
tionship would  then  imply  that  the  decrease  in  peak 
pressure  measured  for  successive  shots  can  be  attributed  to 
the  removal  of  the  gold  and  silicon  dioxide  films  by  suc- 
cessive shots.  Unfortunately,  the  computer  model  did  not 
predict  this  decrease. 
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2.  Pressure  Measurements  at  Atmospheric  Pressure 

As  has  been  previously  described  the  measurements 
of  pressure  pulses  generated  at  atmospheric  pressure  were 
of  small  amplitude.  These  results  were  as  expected  since 
previous  work  by  Magee  and  others  [Ref.  29]  noted  that  in 
comparison  to  pressure  pulses  measured  in  vacuum,  pressure 
pulses  measured  in  atmosphere  on  identical  targets  were 
significantly  smaller. 

In  the  atmosphere  measurements  performed  in  the 
current  investigation  the  pulses  measured  were  generally 
of  less  than  10  bars  magnitude.  In  only  one  case,  the 
1000/2000  system  was  there  any  variation  of  the  pressure 
measured  with  shot  number.  The  pressures  measured  on  the 
thicker  films,  namely  the  2000/3000  system  and  the  0/5000 
system,  remained  almost  constant  throughout  the  data  run 
with  average  values  of  7.4  and  5.4  bars  respectively. 

For  the  data  from  the  1000/2000  system  a linear 
regression  analysis  identical  to  that  performed  on  the  data 
from  the  vacuum  cases  was  performed.  However,  no  corrections 
were  necessary  for  absorption  by  the  entrance  window  since 
none  was  used.  The  linear  regression  analysis  yielded  a 
value  of  k equal  to  -0.046  and  a correlation  to  the  expo- 
nential relation  between  pressure  and  shot  number  of  0.87. 
This  indicates  that  while  an  exponential  relationship 
between  shot  number  and  pressure  probably  exists,  is  is  a 
rather  weak  relationship,  much  weaker  than  that  for  the 
vacuum  case. 
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The  laser-supported  detonation  wave  model  for  the 
generation  of  the  pressure  pulse  leading  to  equation  10 
was  used  to  attempt  to  predict  the  pressure  pulse  using  the 
experimenter-generated  computer  code  previously  discussed. 
The  results  were  uniformly  dismal.  Comparison  of  the  values 
tabulated  in  Tables  VI,  VII,  and  VIII  show  that  the  pre- 
dicted pressures  were  consistently  at  least  two  orders  of 
magnitude  smaller  than  the  experimentally  measured  pulses. 
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VI.  CONCLUSIONS 


A.  DEPENDENCE  OF  PRESSURE  UPON  SHOT  NUMBER  AND  FILM  THICKNESS 
Analysis  of  the  curves  of  pressure  versus  shot  number 
shown  in  Figures  13,  15  and  21  indicates  that  there  exists 
an  exponential  relationship  between  the  pressure  generated 
and  the  shot  number.  Since  total  removal  of  the  silicon 
dioxide  and  gold  films  covering  the  gallium  arsenide 
substrate  of  the  solar  cell  model  was  noted  at  some  point 
during  the  data  runs  it  is  concluded  that  there  exists  some 
relationship  between  the  peak  pressure  generated  and  the 
remaining  thickness  of  the  thin  layers  of  silicon  dioxide 
and  gold  remaining  on  the  gallium  arsenide  substrate.  The 
fact  that  the  gallium  arsenide  substrate  was  entirely 
exposed  at  some  points  after  only  six  shots  on  the  1000/2000 
system  at  atmospheric  pressure  and  that  after  this  exposure 
was  noted  the  pressure  generated  remained  relatively 
constant  seems  to  validate  this  conclusion.  It  is  unfortu- 
nate that  no  similar  observations  could  be  made  during  the 
vacuum  cases  but  the  position  of  the  solar  cell  model  in 
the  test  chamber  and  the  opacity  of  germanium  to  visible 
light  precluded  observation  of  the  damage  to  the  target 
until  completion  of  the  data  run. 
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B.  DETERMINATION  OF  THE  PRESSURE  GENERATION  MECHANISM 
The  primary  methods  used  to  attempt  to  ascertain  the 
mechanism  whereby  the  observed  pressure  pulses  were 
generated  were  examination  of  the  duration  of  the  measured 
pulse  and  a comparison  of  the  measured  pulse  magnitude  with 
that  predicted  by  whichever  of  the  two  models  considered 
in  the  previous  sections  of  this  report  was  applicable  to 
the  conditions  of  the  particular  experiment.  . 

The  comparison  of  the  initial  pressure  pulses  measured 
in  vacuum  to  those  predicted  by  the  plasma  tamping  model  of 
equation  30  yielded  reasonable  agreement  between  the  observed 
and  predicted  pulse  amplitudes.  However,  the  consistent 
underestimation  of  the  pulse  amplitude  by  equation  30  indi- 
cates that  some  other  mechanism  may  possibly  be  contributing 
significantly  to  the  generation  of  the  pressure  pulses. 

In  an  attempt  to  determine  what  this  possible  alternative 
mechanism  might  be  an  examination  of  the  duration  of  the 
measured  pressure  pulse  was  made.  In  general,  the  duration 
' of  the  pulses  remained  near  400  nanoseconds.  This  is 
about  one  and  one-half  the  duration  of  the  250-nanosecond 
pulse  of  the  Carbon  Dioxide  laser  used  to  irradiate  the 
solar  cell  model.  Yang  [Ref.  28]  has  reported  that  the 
pressure  pulse  is  expected  to  have  a duration  close  to  that 
of  the  irradiating  laser  for  pressure  pulses  generated  by 
plasma  tamping.  This  fact  gave  strong  indication  that 
possibly  some  other  mechanism  than  plasma  tamping  might  be 
the  dominant  mechanism. 

■ 
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A further  indication  to  this  effect  is  the  requirement 

that  power  densities  on  the  order  of  10  watts/cm  are 

needed  for  plasma  production  [Ref.  41],  The  power  density 

used  throughout  the  current  investigation  was  of  the  order 
7 2 

of  10  watts/cm  , much  lower  than  the  power  density  needed 
for  plasma  production.  This  strongly  suggested  that  a 
plasma  was  not  being  formed  by  laser  irradiation  of  the 
target,  thereby  strongly  suggesting  that  the  plasma  tamping 
model  was  not  wholly  applicable  to  the  conditions  of  the 
experiment. 

Based  on  these  two  observations  the  plasma  tamping  model 
was  discarded  as  the  chief  mechanism  for  pressure  generation 
in  vacuum  under  the  conditions  of  the  current  investigation. 

Ti.o  applicability  of  the  laser-supported  detonation  wave 
model  was  not  considered  in  the  vacuum  case  since  the  forma- 
tion of  such  a detonation  wave  required  pressures  far  in 
excess  of  the  10  5 torr.  pressure  used  in  the  vacuum  measure- 
ments. This  model  was  examined  as  a possible  mechanism  for 
the  results  obtained  at  atmospheric  pressure  but  again  the 
fact  that  the  power  density  needed  for  plasma  production 
was  not  met  forced  its  exclusion  as  a possible  mechanism. 
Also,  the  uniformly  poor  prediction  of  the  model  in  the 
computer  code  would  have  forced  its  discard  regardless  of 
whether  or  not  the  conditions  for  plasma  production  were  met. 

The  final  mechanism  examined  as  a possible  explanation 
of  the  pressure  pulse  generation  was  the  thermomechanical 
mechanism  described  by  Ready  [Ref.  42],  This  mechanism 
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appeared  attractive  since  it  did  not  require  removal  of 
material  from  the  target  or  the  formation  of  a plasma.  All 
that  is  required  is  rapid  enough  heating  of  an  absorbing 
medium  to  cause  the  thermal  expansion  which  generates  the 
pressure  wave.  The  mechanism  suggests  that  the  duration  of 
the  pressure  wave  generated  is  not  tied  to  the  laser  pulse 
length  but  to  the  time  required  for  the  absorbing  medium  to 
thermally  expand  to  its  greatest  extent  and  than  to  relax 
to  its  ambient  thickness.  In  general,  this  time  would 
probably  be  much  longer  than  the  laser  pulse. 

The  strong  possibility  that  this  mechanism  is  the  dominant 
mechanism  in  the  current  investigation  is  supported  by  two 
factors.  The  first  is  the  long  duration  of  the  pressure 
pulse  when  compared  to  that  of  the  laser  pulse  and  the  second 
is  that  the  pressure  pulse  measured  for  the  very  thick  0.5000 
system  was  of  much  smaller  amplitude  and  greater  duration 
than  those  measured  for  the  thinner  1000/2000  and  1000/1000 
systems.  This  small  amplitude  and  long  duration  is  attributed 
to  the  ability  of  the  thick  film  to  absorb  and  dissipate 
the  thermal  expansion  over  a much  greater  thickness  than 
the  thin  film.  This  accounts  for  the  long  duration  of  the 
pressure  pulse  for  the  thick  film  since  it  takes  a longer 
time  for  the  thick  film  to  expand  than  for  a thin  one. 
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C.  DETERMINATION  OF  POWER  THRESHOLD  FOR  REMOVAL  OF 

COVERING  FILMS 

In  all  cases  except  the  0/5000  system,  the  layers  of 

silicon  dioxide  and  gold  were  found  to  have  been  rapidly 

7 2 

removed  by  the  relatively  low  power  density  of  10  watts/cm 
used  throughout  the  current  investigation.  In  the  cases 
where  observation  was  made  of  the  number  of  shots  necessary 
for  this  complete  removal  to  occur,  the  atmospheric  irradia- 
tions, it  was  observed  that  no  more  than  six  shots  were 
required  before  gallium  arsenide  was  visible  at  some  point 
on  the  surface  of  the  solar  cell  model. 

If  the  dominant  mechanism  at  work  is  the  thermomechanical 
mechanism  described  in  the  previous  section,  this  removal 
of  the  films  of  gold  and  silicon  dioxide  is  attributed  to 
spallation.  A portion  of  the  pressure  wave  generated  by 
thermal  expansion  of  the  gold  layer  is  reflected  at  the  gold- 
gallium  arsenide  interface  and  is  of  sufficient  magnitude 
to  exceed  the  strength  of  the  bond  holding  the  silicon 
dioxide  to  the  gold  thereby  causing  spalling  of  the 
silicon  dioxide  layer.  X similar  sequence  of  reflected 
waves  could  account  for  the  removal  of  the  gold  coating. 

It  should  be  noted  that  this  explanation  for  the  removal 
of  the  covering  films  was  hot  verified  by  experimental 
work  during  the  current  investigation  but  is  merely  offered 
as  a possible  hypothesis  to  explain  the  film  removal. 
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D.  SUMMARY  OF  RESULTS 


The  main  conclusions  arrived  at  in  the  examination  of 
the  results  of  the  current  investigation  are: 

1.  For  films  of  gold  less  than  5000  Angstroms  thick, 
there  is  a relationship  between  the  pressure  generated  and 
the  shot  number  of  the  form 


V I Y 

P = P d • + P 
o 


The  value  of  k*  is  larger  in  vacuum  than  for  irradiation  at 
atmospheric  pressure  and  the  value  of  P in  atmosphere  is 
probably  very  much  smaller  than  the  corresponding  value 
in  vacuum.  Since  removal  of  the  gold  and  silicon  dioxide 
layers  was  observed  after  a series  of  shots,  a dependence 
of  the  constant  k'  upon  the  initial  thicknesses  of  the  films 
is  indicated. 

2.  The  principal  mechanism  for  the  generation  of  the 
pressure  pulses  observed  is  apparently  thermomechanical 
and  due  to  thermal  expansion  of  the  gold  layer. 

3.  The  coatings  of  silicon  dioxide  used  as  an  anti- 
reflective  coating  on  the  solar  cell  model  are  highly 

vulnerable  to  irradiation  by  a Carbon  Dioxide  laser  even 

7 2 

at  the  low  power  density  of  10  watts/cm  used  in  the 
current  investigation.  However,  this  vulnerability  may  be 
significantly  reduced  by  making  the  thickness  of  this  layer 
much  greater.  The  small  amount  of  data  obtained  for  the 


67 


irradiation  of  a 5000  Angstrom  thick  gold  layer  indicates 

that  the  contact  grid  composed  of  a gold  layer  of  this 

thickness  is  able  to  absorb  power  densities  on  the  order  of 
7 2 

10  watts/cm  without  damage. 
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VII . RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

Several  additional  lines  of  investigation  on  the  effects 
of  laser  irradiation  on  solar  cells  are  indicated  by  the 
results  of  the  current  investigation.  The  first  of  these 
areas  is  verification  of  the  results  reported  herein. 

These  conclusions  were  based  on  a very  small  data  base  and 
hence  are  more  indications  of  trends  rather  than  fully 
justified  arguments.  However,  before  this  can  be  done  the 
experimental  apparatus  requires  modification  in  some  way  to 
permit  data  to  be  taken  under  vacuum  conditions  without 
destruction  of  the  entrance  window  to  the  test  chamber. 

Also,  several  other  sets  of  film  thicknesses  should  be 
examined  since  the  current  investigation  shows  that  the  1000 
and  2000  Angstrom  thicknesses  of  silicon  dioxide  are  too 
thin  to  have  a significant  effect  on  the  pressure  generating 
mechanism. 

Another  possible  area  of  investigation  would  be  an 
attempt  to  time -resolve  the  pressure  pulse  generated,  par- 
ticularly to  determine  its  time  of  onset  and  to  attempt  to 
resolve  the  contributions  of  each  layer  coating  the  gallium 
arsenide  to  the  total  pressure  pulse  measured. 

Finally,  since  the  results  of  this  experiment  seem 
to  indicate  thermomechanical  generation  of  the  pressure 
wave,  investigation  should  be  conducted  to  verify  that  this 
is  indeed  the  mechanism  operating.  A refinement  of  the 


computer  model  used  in  the  current  investigation  to  accomplish 
this  is  indicated.  Also  an  attempt  to  predict  the  pressure 
pulse  as  a function  of  time  should  be  attempted  utilizing 
the  material  presented  in  Appendix  C. 


70 


J 


k 


APPENDIX  A 


TABLES 


TABLE  I 


NPS  CC>2  TEA  Laser  Characteristics 


Energy 


-up  to  18  joules* 


Pulse  Width-  -------------  -250  nanoseconds 


Peak  Output  Power- 


up  to  25  megawatts* 


Wavelength 
Beam  Shape 


10.6  ym 


Marx  Input  Voltage 


■square,  4.5  cm  x 4.5  cm 
-------  66  kilovolt 


Laser  Efficiency 


•5%  minimum 


a.  Only  14.4  joules  have  been  produced  with  a peak  power 
of  about  21  Megawatt. 

b.  Considers  only  electron  produced  peak. 

c.  Unfocused,  75  cm  from  output  mirror. 


TABLE  II 


HADRON  Ballistic  Thermopile 

Ballistic  Time  Constant-  ----------  -50  seconds 


Equivalent  Noise  Level 
Maximum  Energy  Input  - 


300  microjoules 
- - 1000  joules 


Maximum  Power  Input - 


10  watts 


Sensitivity- 


44.9  microvolts/ joule 
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W 


Calculated;  see  Appendix 


TABLE  IV 


Results  for  Data  Run  in  Vacuum  for 
a 1000  A Layer  of  Si02  on  a 1000  A of  Au 


Shot 

Peak  Pressure  (bar) 

Pulse  Length 

Number 

Uncorrected 

Corrected 

(nsec. ) 

1 

2 

1440 

2100 

400 

3 

518 

1016 

300 

4 

C 

232 

485 

400 

J 

6 

144 

301 

240 

7 

217 

452 

400 

8 

144 

301 

400 

' 9 

130 

272 

340 

10 

80 

166 

320 

11 

97 

203 

400 

12 

97 

203 

400 

13 

14 

130 

272 

300 

IS 

ns 

239 

280 

16 

100 

209 

320 
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TABLE  V 


Results  for  Data  Run  in  Vacuum  for  a 
1000  A Layer  of  SiO^  on  a 2000  A Layer  of  Au 


Shot 

Peak  Pressure  (bar) 

Pulse  Length 

Number 

Uncorrected 

Corrected 

(nsec . ) 

1 

480 

2 

462 

764 

440 

3 

406 

791 

440 

4 

317 

661 

440 

5 

231 

482 

400 

6 

7 

289 

444 

440 

/ 

8 

259 

540 

320 

9 

144 

301 

240 

10 

115 

239 

280 

74 


TABLE  VI 

Results  for  Data  Run  in  Atmosphere  for 
a 1000  A Layer  of  Si02  on  a 2000  A Layer  of  Au 


Shot 

' 4 

Peak  Pressure 

Pulse  Length 

Number 

(bar) 

(nsec . ) 

1 

2 

14 

500 

3 

13 

600 

4 

5 

6 

7 

14 

500 

8 

14 

500 

9 

12 

-- 

10 

-- 

11 

12 

13 

12 

600 

14 

9 

600 

15 

9 

600 

16 

9 

500 

17 

7 

500 

18 

6 

500 

19 

9 

500 
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TABLE  VII 


Summary  of  Results  for  2000/3000  and  0/5000 
Systems  at  Atmospheric  Pressure 


Total  Number 

Mean  Peak 

Mean  Pulse 

System 

of  Shots 

Pressure 

Length 

(Bar) 

(nsec. ) 

2000/3000 

14 

8.0 

600 

0/ 500a 

7 

5.4 

2000 

a.  Pressure  pulses  for  this  system  were  of  such  small 

amplitude  that  it  was  extremely  difficult  to  separate 
them  from  system  noise.  Therefore,  the  peak  pressure 
and  pulse  duration  upon  which  these  mean  values  are 
based  are,  at  best,  estimates  of  true  values. 
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TABLE  VII 


ects  change  in  reflectance  when  entire  silicon  dioxide  layer  has  been  remove 


TABLE  IX 

Peak  Pressure  for  1000/1000  System  in  Vacuum  Corrected 
for  an  Average  Saturation  Pressure  of  215  Bar. 


Shot 

Saturated  Peak 

Peak 

Pressure  Corrected 

Number 

Pressure  (Bar) 

for 

Saturation  (Bar) 

1 

2 

2100 

1885 

3 

1016 

801 

4 

r 

485 

270 

5 

6 

301 

86 

7 

452 

237 

8 

301 

86 

9 

272 

57 
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TABLE  X 


Peak  Pressure  for  1000/2000  System  in  Vacuum  Corrected 
for  an  Average  Saturation  Pressure  of  215  Bar 


Saturated  Peak 


Peak  Pressure  Corrected 


Pressure  (Bar) 


for  Saturation  (bar) 


APPENDIX  B 

Calculation  of  Values  of  nR,  nj  and  k for  Gallium 
Arsenide,  Gold  and  Silicon  Dioxide 

From  Ref.  43  the  values  of  the  real  and  imaginary 
parts  of  the  index  of  refraction  are  given  by 


„ _ fY+ea/2 

nR  ’ ^ 


„ - fY-0>l/2 

nI  " (-r“) 


where 


6 = l-Cu>pV+*c2)) 


y = [B2+(1-3)2(vc/u))2]1/2 


in  which  w is  the  laser  angular  frequency,  u>  is  the  plasma 
frequency  given  by 


cop  =[4Tre2Ne/m]1/2 

= 8.5  (pz)1/2*1015  (radians/sec) 


where  p is  the  density  and  z is  the  number  of  free  electrons 


per  atom. 
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The  quantity  vc  is  the  collision  frequency  and  is 


Vc  - up/4ir«0 


in  which  crQ  is  the  static  conductivity. 

The  absorption  coefficient,  k,  is  defined  by 


k = (2a)/c)nj 


where  c is  the  speed  of  light  in  vacuum  (S.O'IO1  cm/sec). 
The  "skin  depth"  is  defined  by 


6 = k'1  = XQ/4TTnI  . 


(B-8) 


For  the  CO.,  laser,  X = 10.6  urn  and  cu  = l.yS'lO"''4  radians/ 
2 o 

sec  ^ . 

The  transmittance,  the  ratio  of  transmitted  to 
incident  energy,  is  defined  by 


T = e 


•kx 


(B-9) 


where  x is  the  thickness  of  the  sample. 
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1.  Calculation  of  n^,  n^  and  k for  gold. 

For  gold  z = 1 and  aQ  = 3.7*10 17  sec  1 [Ref.  37] 


to  = 8.4  (pz)1/2-1015  sec'1 
P 

= S^-IO16  radians/sec  1 


v = uj  /4Tra 
c p7  o 

^ iri14  -1 
= 2.9*10  sec 


Then,  substituting  and  vc  into  equations  B-3  and  B-4, 
it  follows 


6 = 1-  (cOp//co2  + v2))2 


= -1.2*10 


Y = [B2  + (1-3) 2 (vc/6)2]1/2 


= 2.3*10  . 


Substituting  these  quantities  into  equations  B-l  and  B-2, 
the  values  for  the  real  and  imaginary  parts  of  the  index 


of  refraction  are 
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nR  = 74.2 
n j = 132.3 

From  equation  B-7,  the  value  of  the  absorption  coefficient 
for  gold  is  calculated  to  be  k = 1.5*10^  cm 
2.  Calculation  for  nj  and  k for  silicon  dioxide. 

Reference  53  shows  the  transmittance  of  silicon  dioxide 
v-  at  10.6*1,111  to  be  0.82  for  a sample  6680  A thick.  From 
equation  B-9,  it  is  known  that 


T = 


Then 


and 


In  T = -kx 


k = -In  T/x  . 
From  the  parameters  given  in  Ref.  53, 

k = 3.0-103  cm"1 
From  equation  B-7,  nj  is  found  to  be 

nj  = 0.25  . 


(B-10) 
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3..  Calculation  of  Absorption  Coefficient  for  Gallium 
Arsenide . 

From  equation  B-7,  it  is  known  that 


k = (2<jo/c)nj  . 


Table  III  of  Appendix  A shows  that  at  10.6  pm  n^  = 8.8*10 

- 2 - 1 

Therefore  k * 1.0*10  cm  . Note  that  this  gives  a skin 


depth  of  about  1.0  meter  indicating  that  gallium  arsenide  is 
essentially  transparent  to  the  radiation  from  a C02  laser. 


'■ 
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APPENDIX  C 

Development  of  Stress  (cr(x,t))  and  Temperature  (T(x,t)) 
as  Functions  of  Position,  x,  and  time,  t 

An  attempt  was  made  to  calculate  the  stress  generated 
in  the  solar  cell  model  as  a function  of  the  distance  the 
stress  has  propagated  into  the  model,  x,  and  the  time 
elapsed,  t. 

Michaels  [Ref.  39]  has  solved  for  the  stress  generated 
in  a bar  as  a function  of  x and  t. 

For  an  elemental  slab  of  a bar  of  cross-sectional  area 
A,  the  net  force  acting  on  the  slab  is 

A3a/3x  dx  = pA  dx3^u/3t^  (C-l) 


in  one  dimension  where 

u = longitudinal  displacement  of  slab, 
i t = time, 

x * longitudinal  coordinate  of  bar, 
p = mass  density  of  bar, 

a * longitudinal  stress  generated  in  bar. 

In  order  for  a one-dimensional  treatment  of  the  problem 
to  be  justified  the  effects  of  lateral  inertia  and  shear 
must  be  negligible.  The  laser  spot  size  used  in  the 
current  investigation  irradiated  the  entire  surface  of  the 
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target  so  that  the  effects  of  lateral  inertia  and  shear 
can  be  neglected. 


From  thermo-elasticity 


a = E(e  - aT(x,t)) 


(C-  2) 


= elastic  modulus  in  dynes/cm  , 

O 

= coefficient  of  expansion  in  cm/ cm  K, 
= temperature  of  bar, 

= total  longitudinal  strain. 


e = 9u/9x 


and  c = E/p  where  c is  the  stress  wave  velocity,  equation 
C-l  becomes 


92u/9x2  - a9T(x,t)/9x  = l/c292u/9t2  . (C-3) 


For  a finite  bar  extending  from  0 to  L,  we  have  the  boundary 
conditions  that  the  stress  at  each  end  of  the  bar  is  zero. 
That  is 


o|x=0  a E[9u/9x  - aT(x,t)]x=0  = 0 


(C-4) 
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F 


or 


3u/3x|x=q  = aT(0,t)  (C-5) 

311/  Ox  | v_l  = aT(L,t)  (C-6) 

The  initial  conditions  for  the  bar  are 


u|t=0  = 0 and  3u/3t|t=0  = 0 

Recall  equation  C-3.  Taking  its  finite  Fourier  cosine  transform 


l/c232/3t2  j u(x,t)  cos  (mrx/L)  dx 
o 

L 

* j 32u/3x2  cos  (mrx/L)  dx 


a J 3T(x,t)/3x  cos  (mtx/L) 


(C-7) 


Integrating  the  first  term  on  the  right-hand  side  of  this 
equation  by  parts,  the  result  is 


jj 

j 32u/3x2  cos  (Bpi)  dx  = (-l)n[3u/3x]x=L  - [3u/3x]x=( 


- n 2 u(x,t)  cos  (^p-)  dx 


Invoking  the  boundary  conditions,  equations  C-S  and  C-6, 
this  becomes 
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J 32u/3x2  cos  (SJi)  dx  = C-l)n  ctT(L,t)  - aT(0,t) 


2 2 L 

‘ ^~T~  f uC*.t) 


cos  (2H)  dx 


(C 


Using  equations  C-5,  C-6  and  C-8  equation  C-7  becomes 


1/c2  32uc(n,t)/3t2  = (-l)naT(L,t)  - aT(0,t) 


- n2ir2/L2uc (n,t)  - a9  (n,t) 


(C 


where 


uc(n,t)  = J u(x,t)  cos  (nirx/L) 


dx 


CC 


and 


>c(n,t)  = J~  3T(x,t)/3x  cos  (nirx/LJ 


dx 


(C 


Equation  C-9  reduces  to 


32uc(n,t)/3tZ  + (nZTTZcZ/L“)uc(n,t)  = cZa  [ (- 1)  nT  (L , t)  - T(0, 


2_2_2 


- 0c(n,t)] 


CC 


Now,  define 


2 2-  2 2 /t2  „ , 
a)  = n /I  c /L  and 


P(n,t)  - c2ct[(-l)nT(L,t)  - T (0 , t)  - 0c(n,t)] 


-8) 

-9) 

-10) 

-ID 

t) 

-12) 
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and  equation  C-12  becomes 


32uc/3t2  + oo2uc  = P(n,t)  (C-13) 

I 

The  solution  to  this  differential  equation  is  [Ref.  39] 

■ t 

uc(n,t)  = 1/w  j P(n,z)  sin  (uj(t-z))  dz 
o 

t 

= cLa/nir  /[(- l)nT(L,z)  - T(0,z)  - ?r(n,z)] 
y c 

0 

sin  (mrc/L  (t-z))  dz  (C-14) 

By  integrating  equation  C-ll  by  parts  it  is  found  that 

L 

^(njt)  = T(L,t)(-l)n  - TC0,t)  + mr/L  j~  T[x,t)  sin  (n  x/L)  dx 

o 

Substituting  this  into  equation  C-14  it  is  seen  that 
t L 

uc(n,t)  = -ca  ysin(mrc/L(t-z))  j T(x,z)  sin(mrx/L)  dxdz  (C-15) 
o o 

By  using  the  definition  of  the  finite  Fourier  transform  on 
equation  C-15  it  is  found  that 
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u(x,t) 


uc(0,t) 


l l cosSJS  co,  f 3in2IH^ii42  /TCx,z)sin2Ji  dx 


u„(0,t)  = cct  j sin(O)  J T(x,z)  sin(O)  dxdz  = 0 


u(x,t)  .l2“zcos(!l”)  fT(x,z)sln(2fi)dxd; 

L n=l  h J L J L 


CC-16) 


Hence,  if  T(x,t)  is  known,  the  stress  in  the  bar  can  be 
obtained  from  equation  C-2.  That  is 


a(x,t)  = pc  [3u(x,t)/3x  - aT(x,t)] 


(C-17) 


The  problem  of  determining  T(x,t)  for  a two-layered 
wall  analogous  to  the  silicon  dioxide  and  gold  layers  of 
the  solar  cell  model  used  in  the  current  investigation  has 
been  solved  by  Griffith  and  Horton  [Ref.  40]. 
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The  solution  assumed  that  the  area  of  the  wall  was 
sufficiently  large  that  the  heat  is  conducted  perpendicularly 
to  its  surface,  which  is  taken  as  the  plane  x = 0.  The 
wall  consists  of  two  components,  one  of  thickness  a,  extending 
from  the  plane  x = 0 to  x = a,  while  the  second  component 
extends  from  x = a to  x = 

The  temperature  T(x,t)  must  satisfy  the  following 
boundary  value  problem 

3T/3t  = D32T/3x2  t > 0,  a > x >_  0;  D = Dx  (018) 

t > 0,  °°>^x>_a;  D = D 2 


T1(x,0)  = 0 

(019) 

T1(a,t)  = T2(a,t) 

(C- 20) 

kl3Tl/3xlx-a  * VV3xi 

^ x=a 

(C- 21) 

(-k13Tl/3x)|J(.0  = H for 

o 

A 

r+ 

A 

8 

(C- 22) 

T2(«>,t)  = 0 for  0 < t < 

00 

( C - 2 3 ) 

where  subscript  1 refers  to  the  top  layer  1,  subscript  2 
refers  to  the  base  layers  2 and 
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x = distance  perpendicular  to  face  of  wall, 
t = time, 

D = thermal  diffusivity  = k/S, 
k = thermal  conductivity , 

H = constant  flux  of  heat  perpendicular  to  the  wall 
surface , 

S = volume  specific  heat  (specific  heat  times  density). 
The  solutions  to  this  boundary  value  problem  are  found  to 
be : 

Ti(x,t)  = |L  e-x2/4Dit  _ xU.er£{x/2^T])] 

- I ”(-l)n[5fgTE'(-c-{x>2a(n>l)}2/4D1t  + {x- 2a  (n+1)  ^MD^) 

Y n=l  Y ” 

-{x+2a(n+ l)Hl-erf(x+2ct(n+1))> 

2/d^T 

+ {x-2ct  (n+1)  Hl-erf  [-  (~c*2a  ) J } ] } (C-23) 

2/D^T 


2HXK  “ (-l)n{2|°2^,c-{x-a{l- JP2/D1(2n-<-l)}}2/4D2t 


T2(x,t)  y n=Qv  y'  rr 


’■)  x-a{l- ^d7/d7 (2n+l)  } 

-[x-a{l  -4/  (2n+l)  }]  (1-erf  [ ])} 

U1  2Vd~T 


])}  (C- 24) 


The  surface  temperature  TQ(t)  is  given  by  setting  x = 0 
in  equation  C-24.  Then 


ToCt)  = H/k1[2^£'-  is  (^)n{GTe‘a2fn+1)2/tDl 


a(n+l)(l-er£[^il])}] 

Vd 7T 


(C-25) 


As  can  be  seen,  by  applying  equations  C-24  and  C-25 
to  equation  C-17  it  is,  in  principle,  possible  to  predict 
the  stress  in  the  solar  cell  model  as  a function  of  both 
position  and  time.  A computer  code  to  accomplish  this  is 
available  from  Dr.  Tom  Magee  of  Stanford  Research  Institute. 
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APPENDIX  D 

NPS  C02  Laser  Operating  Procedure 

The  following  recommended  procedure  includes  changes 
made  to  the  operating  procedure  given  by  Strickland 
[Ref.  30]  and  is  included  as  a reference  for  future  users 
of  the  laser  system.  The  procedure  is  as  follows: 

1.  Open  main  valves  at  gas  tanks  and  set  flow  rate 
at  following  values: 

He  = 4.0 
C02  * 5.0 
N2  = 4.0 

Allow  a flow  of  four  minutes  before  firing  the  laser. 

2.  Set  trigger  generator  to  warm-up  position  and  allow 
a two-minute  warm-up  period.  Then  set  trigger  generator 

to  single  shot  position. 

3.  Turn  interlock  switch  on  Circuit  8 to  the  "ON" 
position.  This  activates  the  red  warning  lights  at  entrances 
to  the  laboratory  and  puts  power  to  the  DC  power  supply  of 
the  laser. 

4.  Ensure  that  the  north  door  to  the  laboratory  is 
locked  to  prevent  unwary  individuals  from  entering  by  this 
means . 

5.  Turn  on  DC  power  supply  utilizing  the  switch  located 
on  front  cover.  This  will  light  both  the  amber  and  green 
lights  on  the  power  supply.  Allow  a two-minute  warm-up 
period. 
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6.  Put  on  safety  goggles  with  plexiglass  window. 

7.  Set  Marx  generator  spark  gap  pressure  corresponding 
to  desired  firing  voltage  following  posted  chart. 

8.  Set  desired  Marx  generator  discharge  voltage  on 
variac  at  the  DC  power  supply. 

9.  Set  Marx  generator  charging  voltage  on  preset 
voltage  meter  at  the  DC  power  supply.  The  voltage  set 
should  be  one-third  of  the  Marx  discharge  voltage. 

10.  Initiate  charging  of  the  system  by  pressing  the 
remote  charge  switch  mounted  next  to  the  trigger  generator 
or  by  pressing  the  "CHARGE"  switch  at  the  DC  power  supply. 

This  will  extinguish  the  green  light  and  activate  the  red 
light  on  the  power  supply. 

11.  When  the  preset  voltage  level  is  reached,  the  audible 
alarm  will  be  activated  and  the  red  light  will  be  extinguished. 
Firing  of  the  laser  can  now  be  accomplished  by  pushing  the 
button  on  the  trigger  generator.  If  the  laser  is  not  fired 
within  five  seconds  after  completion  of  charging  the  auto- 
matic dump  operates  and  the  system  is  discharged  to  ground. 

12.  If  necessary,  the  system  can  be  dumped  at  any  time 
after  initiation  of  charging  either  with  the  remote  dump 
switch  at  the  trigger  generator  or  the  local  dump  switch 
at  the  DC  power  supply. 

13.  Further  firing  at  the  preset  voltage  can  be 
accomplished  by  repeating  the  procedure  beginning  at  step  8 
or  at  another  voltage  beginning  at  stop  5. 


1 
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NPS  CO^  Laser  Alignment  Procedure 

The  following  alignment  procedure  is  suggested  for 
future  users  of  the  C02  Laser  System.  It  should  be  imple- 
mented whenever  the  output  of  the  laser  is  found  to  have 
significantly  changed  or  when  the  Sodium  Chloride  (NaCl) 
Brewster  windows  are  installed  after  polishing.  The 
alignment  procedure  requires,  in  addition  to  the  C02 
laser  itself,  the  following  equipment: 

1.  Helium-Neon  (HeNe)  CW  laser  with  laser  aligner 
attachment. 

2.  C02  laser  alignment  blocks  (2) 

3.  HADRON  Ballistic  Thermopile 

4.  Millimicrovoltmeter 

5.  Copper  Mirror  (f  = 1.0  m) . 

The  alignment  procedure  assumes  that  the  NaCl  windows  are 
not  installed  at  initiation  of  the  procedure.  It  is  broken 
into  two  sections:  rough  alignment  and  fine  tuning. 

Referring  to  Figure  25,  the  procedure  is  as  follows: 

A.  Rough  Alignment 

1.  Remove  Germanium  exit  mirror  and  its  mount. 

2.  Install  laser  alignment  blocks. 

3.  Install  NaCl  windows. 

4.  Mount  HeNe  laser  about  1.0  meter  in  front  of  output 
end  of  C02  laser. 
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5.  Adjust  beam  from  HeNe  laser  so  that  it  passes 
through  the  holes  in  center  of  laser  alignment  blocks. 

6.  Place  a sheet  of  Mylar  between  output  end  of  CO 2 
laser  and  HeNe  laser. 

7.  Adjust  Copper  Mirror  (f  - 3.0  m)  at  back  end  of 
CO^  laser  resonant  cavity  until  the  beam  from  the  HeNe 
laser  reflected  by  this  mirror  passes  back  through  holes 
in  alignment  blocks  onto  the  Mylar  sheet. 

8.  Adjust  the  Copper  Mirror  until  the  reflected  beam 
is  coincident  with  the  beam  from  the  HeNe  laser  incident 
on  the  Mylar. 

9.  Reinstall  Germanium  output  mirror  and  mount. 

10.  Adjust  germanium  mirror  until  the  HeNe  beam  reflected 
from  this  window  onto  the  Mylar  sheet  is  coincident  with 

the  incident  beam. 

11.  The  CO^  laser  is  now  rough  aligned.  Remove  the 
laser  alignment  blocks.  This  will  require  removal  of  the 
NaCl  windows.  When  removing  these  windows  do  not  change 
their  orientation. 

12.  Reinstall  NaCl  windows. 

13.  Remove  HeNe  laser  and  Mylar  sheet. 

14.  Place  a sheet  of  laser  "footprint"  paper  about 
20  cm  in  front  of  the  germanium  output  window. 

15.  Light  off  the  C02  laser  system. 

16.  Starting  at  about  30  kV  discharge  voltage  from  Marx 
generator  fire  the  laser.  If  no  beam  footprint  is  observed 


on  the  footprint  paper,  increase  the  discharge  voltage  in 
increments  of  3 kV  until  a footprint  is  seen.  Do  not 
exceed  45  kV  at  this  stage  of  alignment.  When  footprint  of 
laser  beam  is  seen  proceed  to  step  17. 

17.  Remove  footprint  paper  and  begin  fine  tuning  the 
laser. 

B.  Fine  Tuning 

18.  Install  1.0  meter  focal  length  copper  mirror  about 
1.0  meter  from  output  window  of  laser. 

19.  Connect  input  lead  of  voltmeter  to  calibration 
terminals  of  ballistic  thermopile. 

20.  Mount  ballistic  thermopile  about  1.0  m from  copper 
mirror  of  step  18. 

21.  Adjust  copper  mirror  so  that  beam  from  CC^  laser 
enters  INPUT  cone  of  ballistic  thermopile. 

22.  Using  the  voltmeter  to  monitor  the  total  laser 
energy  output  adjust  the  mirror  at  back  end  of  laser 
resonant  cavity  until  laser  energy  output  is  at  a maximum. 

A wait  of  about  two  minutes  between  laser  firings  will 
yield  best  results. 

23.  Repeat  step  22  with  germanium  output  mirror  of  C02 
laser.  This  completes  the  alignment  procedure.  A calibra- 
tion curve  for  laser  output  as  a function  of  discharge  voltage 
may  now  be  obtained. 
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APPENDIX  F 

COMPUTER  CODE  TO  CALCULATE  PEAK  PRESSURE  PULSES 

THIS  PROGRAM  IS  DESIGNED  TO  CALCULATE  PEAK  PRESSURE 
PULSES  GENERATED  IN  A GALLIUM  ARSENIDE  SCLAR  CELL  MODEL  AS 
THE  SILICON  DIOXIDE  AND  GOLD  LAYERS  ARE  SUCCESSIVELY  RE- 
MOVEC  IN  50  ANGSTROM  DECREMENTS.  THE  CCOE  FIRST  CALCULATES 
THE  FEFLECTANCE  ANO  ABSORPTANCE  OF  THE  MODEL  FOLLOWING  THE 
METHOD  GIVEN  BY  HEAVENS  (REF. 36).  THE  AeSORPTANCE  SO- 
CALCLLATED  IS  THEN  APPLIEO  TO  THE  LASER  ENERGY  OUTPUT  (EN, 
END  TO  GET  THE  ENERGY  OENSITY  INCIDENT  ON  THE  TARGET  (EC. 
£01).  THE  VALUE  CF  ENERGY  DENSITY  IS  THEN  APPLIED  TO  GET  THE 
PRESSURE^ GENE RAT ED  IN  BOTH  ATMOSPHERIC  AND  VACUUM  PRESSURE 

CCNCJhECV?RIOUS  LAYERS  OF  THE  SOLAR  CELL  MODEL  ARE 
IDENTIFIED  BY  A NUMERICAL  DESIGNATION  AS  FOLLOWS: 


1.  AIR 

2.  SILICON  DIOXIDE 

3.  GCLO 

A.  GALLIUM  ARSENIDE 

5.  CUARTZ  GAUGE  (SILICON  CIQXIOE). 

THE  INITIAL  THICKNESS  (D)  AND  THE  REAL  AND  IMAGINARY 
PARTS  OF  THE  INOEX  OF  REFRACTION  (RR  ANC  Ri ) FOR  EACH  LAYER 
ARE  READ  IN  WITH  A MATRIX  NOTATION.  ANOTHER  MATRIX  IS  USED 
TO  ENTER  THE  VALUES  OF  DENSITY  (RHO)  ANC  SOUND  SPEEC  (V)  FOR 
EACH  MATERIAL.  TWO  MATRICES  ARE  USED  SINCE  THE  VALUES  OF  RP 
ANC  PI  ARE  WAVELENGTH  DEPENDENT  ANC  CAN  THUS  BE  CHANGED 
WITHOUT  AFFECTING  THE  VALUES  OF  RHO  ANC  V. 

THE  COMPUTER  PROGRAM  FOLLOWS. 

ENTEF  JCECARCS. 

//CLICKRUN 

JFJ001  FORTGO  ( 2720 ,0420, WS44) » * JACOBS 0 NS MC 1459* 

CIMENSICN  MATRICES  UScO. 

DIMENSION  RR(5)tRI(5)»D(5)»RHC(5)»V(5) 

ENTER  VALUE  OF  WAVELENGTH  IN  METERS.  (WA). 

WA-1C.6E-06 

ENTER  LASER  SPOT  RAOIUS  IN  METERS  (RA). 

R A-6. 5E-03 

ENTER  PLLSE  LENGTH  (TP)  AND  PRESSURE  PULSE  DURATION  AT 
ATMCSFHERIC  PRESSURE  (T).  BOTH  OF  THESE  CUANTITIES  ARE  IN 
UNITS  CF  SECONDS. 

T F-250. 0E-09 
T-5.0E-C7 

ENTER  LASER  ENERGY  OUTPUT  (ENL)  IN  JOULES. 

EM-5.25 

REAC  IN  AND  PRINT  OUT  THE  VALUES  OF  0(1),  RR(I),  ANC  RI(I). 
C(I  ) IS  IN  UNITS  OF  METERS. 

FI-5.14159 

REAC  (5.500)  (RR < I ) , RI ( I ) ,0(1) ,1*1,5) 

5CC  FCFMAT  (3E10.3) 

WRITE  (6,600)  (RR(I),RI(I),D(I>,  1-1,5) 

6CC  FORMAT  (//.10X, 3510.3 ) 

WRITE  (6,601) 

601  FCFMAT  (//) 

REAC  IN  AND  PRINT  OUT  THE  DENSITY  AND  VELOCITY  MATRICES. 


REAC  (5,501)  (RHQ(J).V(J) ,J-2,5) 
FCFMAT  (2E10.3) 

WRITE  (6,604)  ( RHO (K ) , V (K ) , K-2 , 5 ) 


6C4  FORMAT  ( // , 5X , 2E 10. 2) 

WRITE  (6,601) 

CALCULATE  ANC  PRINT  OUT  THE  ACOUSTIC  IMPEDANCES  PCR  USE  IN 
CALCULATING  THE  PRESSURE  IN  VACUUM. 

Z2-RHC ( 2) *V(2 ) 

Z2-PH0 ( 3 ) *V ( 3) 

Z4- RHO ( 4 ) *V (4 ) 
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WRITE  (6*610)  Z2.Z3.Z4 


610  FORMAT 
WRITE 


*601) 


3E10, 


CALCULATE  THE  ENERGY  DENSITIES  IN  JOULES  PER  SQU 
THE  ENERGY  DENSITY  ED  INCLUDES  THE  CORRECTION  FC 
GERMANILM  ENTRANCE  WINDOW  ANC  IS  USEO  FOR  VACUUM 
EOl  CNL  Y CORRECTS  FOR  THE  GERMANIUM  LENS  AND  IS 
ATMOSPHERIC  PRESSURE  CALCULATIONS. 

EN  »0.6$3*ENL 

EM«0.9*ENL 

EC»EN/( (PI )*(RA**2) ) 

ECl*ENl/( (PI )*(RA**2) I 

SP  j 1:14!  I EDl 

WRITE  (6*608)  ED 

606  FORMAT  (5X* 'ENERGY  DENSITY  IS',1E10.3I 
615  FORMAT  (5X.» ATMOSPHERE  ENERGY  1S',1E10.3) 

AT  THISPFcfN|6fHeREFLECTANCE  OF  THE  SYSTEM  OF  T 
IS  BEGUN  TO  BE  CALCULATED.  FIRST  THE  TERMS  INOEP 
FILM  THICKNESS  ARE  CALCULATED  AND  PRINTEC. 


i FOR  THE  GERMANIUM 
SURE  CALCULATIONS. 


SQUARE  METER. 

FOR  THE 
UUM  PRESSURES, 
IS  USEC  FOR 


IS  FCINT  The  reflectance  of  the  system  of  thin  films 

GUN  TO  BE  CALCULATED,  FIRST  THE  TERMS  INDEPENDENT  OF 
THICKNESS  ARE  CALCULATED  AND  PRINTEC. 

A2»(RR( 1)**2)+(RI (i)**2)-(RR(2)**2)-(RI(2)**2) 
E2»((RR (l)*RR(2) )**2)+( (RI ( 1 )+R 1(2)  )**2  ) 
C2«2.0*((RR(ir*RI(2))-(RR(2)*RI(l))) 

A2«(RR(2)**2)  + (P. I (2)**2)-(RR(3)**2)-(RI  ( 2 )**2 ) 
E2*((RR(2)+RR(3)  )**2)  + ((RI(2)+RH2)  )**2) 


,RR(2I*RI (3) )-<RR(3)*RI<2) ) ) 


)**2 ) 


A4«(RR(2)**2)«-(RI(2)**2)-(PR(4)**2)-(RI (4)**2) 
BA«((RR(3) +RR (4) ) **2 ) ♦ ( (RR ( 2) +R  I (4 ) )**2 ) 
C4«2.Q*((RR(3)*RI(4) )-(RR (4)*R I < 2 ) ) ) 

G2-A2/B2 
G2»A2/B2 
GA«A4/B4 
H2*C2/B  2 
H2*C3/B3 
H A*C4/B  A 
WRITE  (6*606) 

6C6  FORMAT  (14X , *G1 ' , SX,  • C-2 ' * 9X , • G3  • ,9X , • G4  • ) 

WRITE  (6.605)  G1*G2*G3*G4 

605  FCRMAT  (10X.4E10.3) 

• WRITE  (6*601) 

WRITE  (6,607) 

607  FCRMAT  ( 14X , »H1' , 9X , • H2* ,9X , • H3 • ,SX , • H4 • ) 

WRITE  (6*605)  Hl*H2,H3,H4 
WFITE  (6*601) 

THE  RACIAL  WAVELENGTH  IS  NOW  CALCULATEC. 

W«(2.0*FI l/WA 

AT  THIS  POINT  THE  TERMS  OF  THE  REFLECTANCE  DEPENDENT  UPON 
THE  THICKNESSES  CF  THE  FILMS  ARE  CALCULATED  ANC  THE  PRESSURE 
FOR  THE  VACUUM  AND  ATMOSPHERIC  CASES  ARE  THEN  COMPUTED. 

THIS  IS  CONE  EY  THE  USE  OF  TWO  DO  LOGPS.  THE  FIRST  STARTS 
WITH  THE  INITIAL  THICKNESS  OF  ThE  SILICON  OIOXICE  LAYER  ANC 
REMOVES  IT  IN  50  ANGSTROM  DECREMENTS.  WHEN  THIS  LAYER  HAS 
BEEN  TOTALLY  REMCVED  THE  COOE  SHIFTS  TO  THE  SECCNC  DC  LOOP 
ANC  REPEATS  THE  PROCEDURE  FOR  THE  GOLC  LAYER. 


CALCULATE  THE  TERMS  OF  THE  REFLECTANCE 
CIOMCE/GCLO  SYSTEM  THAT  ARE  THICKNESS 
DC  15  K-1,50 


CC  15  K-1.50 

E 2*EXP ( ALF2 ) 
C2*C0S( GAM2 ) 
D2*SIN( GAM2) 
P 2»E2*C  2 


F THE  SILICON 
EPENDENT. 


C2-E24C3 

uI’flS!J81t- 


|H3*D3)  )/E2 


Li«((H2AD2)-(G3*03))/E2 
CALCLLAlE  THE  REFLECTANCE,  R,  OF  THE  SYSTEM. 


I 


F23*P3+ (G2*T3)-(H2*U3 ) 
C22*Q3+  H2*T3)+(G2*U3) 
T23«T3+(G2*P3)-(H2*C3) 


R«((T23R*2I+(U23**2))/((P23**2)+(C23**2)  ) 
CALCULATE  THE  ABSCRPTANCE  OF  THE  SYSTEM*  A, 
A«1.0-R 

CACLLATE  THE  PRESSURE  GENERATED  I K THE  GOLD  LAYl 
SOLAP  CELL  MCCEL  FOR  THE  VACUUM  (PBAR)  AND  THE  , 


( PPEAR ) CASES.  UNITS  ARE  I 
PR2*Z3*A*ED/(3.0*TP) 
PF*SCPT (PR2) 


IN  BARS. 


.D  LAYER  OF  THE 
5 THE  ATMOSPHERIC 


PEAR*PR/1.CE05 

PF«(A*EC1)/(V(3)*T) 


PF*  (A*EC1 ) / ( V(  3 . 

FFEAR*PP/1.0E05 

CORRECT  THESE  PRESSURES  FOR  THE  VALUES  EXPECTED  IN  THE 
GALLIUM  ARSENIDE  LAYER  ( PGAAS  AND  PPGAAS  RESPECTIVELY.). 
PGAAS«0.43*PBAR 
FF6AAS*0.43*PPBAR 

PRINT  THE  VALLES  OF  R,  A»  PBAR,  PGAAS  »PPBAR,  ANC  FPGAAS 
CORPESFCNCING  TO  LAYER  THICKNESS  C(2). 

WRITE  ( 6*602 ) D< 2 I ,R» A, PBAR , PGAAS ,PPBAR , PPGAAS 
602  FCFMAT  (5X,7E12.4) 

RECUCE  THE  THICKNESS  OF  THE  SILICON  DICXICE  LAYER  BY  50 
ANGSTROMS  AND  DETERMINE  IF  IT  HAS  BEEN  COMPLETELY  REMOVEC. 
IF  It  HAS  SHIFT  TO  THE  NEXT  CC  LGCP.  IF  NOT  REPEAT  ABOVE 
COMPLTAT ICNS  FOR  THE  NEW  VALUE  OF  0(2). 

C (2  )*C ( 2 )-5.0E~09 
IF  (C(2  i.LE.0.00 ) GO  TO  21 
15  CONTINUE 

21  CC  22  1*1*50 

ALF3*W*C(3)*RI(3) 

GAM3»W*C(3)*RR(3) 

CC3»C0S(GAM3) 

CS3*SIN (GAM3 ) 

E2*EXP ( ALF3 ) 

P4*E3*CC3 


CORRECT  THESE  PRESSURES  FOR  THE  VALUES  EXPECTED  IN  THE 
GALLIUM  ARSENIDE  LAYER  (PGAAS  ANO  PPGAAS  RESPECT! VELY.  ) • 


i PPEAR*  ANC  PPGAAS 


QA*E3*CS3 

T A* ( ( G4*0C3 )+( H4*0S3 ) )/E3 
L4«(<H4*OC3)-(G4*OS3) )/E3 
P34*P4*(G3*T4)-(H3*U4) 
C34*Q4+ (H3*T4)+(G3*U4) 
T34*T4+ (G3*P4)-(H3*C4) 
L34*U4* (H3*P4)+(G3*C4) 
COMPUTE  THE  REFLECTANCE  OF  THIS 


COMPUTE  THE  REFLECTANCE  OF  THIS  NEW  SYSTEM  WHEREIN  THE 
SILICON  CICXICE  LAYER  HAS  BEEN  ENTIRELY  REMOVED. 

RP«((T34**2)+(U34**2) )/((P34**2)4(C24**2)  ) 

CALCULATE  THE  ABSORPTANCE  OF  THE  SYSTEM*  AP. 

AF*1.0-RP 

CALCULATE  THE  VACUUM  ( PP8AR ) AND  ATMOSPHERE  (REAR) 
PRESSURES*  CORRECT  THEM  TO  THE  VALUES  EXPECTEO  IN  THE 
GALLIUM  ARSENIDE  LAYER  (PG  AND  PPGR)  ANC  PRINT  THE  VALUES 
OF  RF*  AF,  PREAR,  PG,  RBAR,  ANC  PPGR  FCR  CORRESPONDING 
THICKNESS  OF  THE  GOLD  LAYER  0(3). 

PR2F«Z4*AP*ED/(3.0*fP) 

PFF*SGRT ( PR2P ) 

PRE AR*PRP/ 1.0E05 
PG*0.43*PRBAR 
PFR*(EDl*AP)/( V(4)*T) 

REAR-PPR/1.0E05 
PFGR*0. 43* RBAR 


TE  (6*602)  D(3)»RP»AP»PRBAR, PG*PBAR»PPGR 
THE  GOLD  LAYER  BY  50  ANGSTROMS  AND  REPEAT  THE 
CNS.  IF  LAYER  HAS  BEEN  COMPLETELY  REMOVEC 


PFGR*0. 4 3* RBAR 
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